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INTRODUCTION 


Jusqu’à présent nous n’avons considéré que des systèmes dont l’énergie peut, en principe, 
augmenter indéfiniment et dont le nombre de micro-états augmente fortement avec l’énergie. 
C’est le cas, par exemple, du gaz parfait classique. De tels systèmes sont appelés normaux et 
la température est une quantité positive. Il existe néanmoins, dans la nature, des systèmes dont 
l’énergie est bornée à la fois inférieurement et supérieurement et pour lesquels la 

quantité ( dE ) peut être négative , dans ce cas, la température absolue devient négative. 

KologLl/ y m 

C’est le cas des systèmes de spin où la projection du spin ne peut prendre qu’un nombre fini 
de valeurs. Pour ces derniers, la température peut être positive ou négative. Si l’on revient au 
cas d’une assemblée de particules n’ayant que deux niveaux d’énergie (associées à deux 
projections du spin), une température négative correspond à ce que l’on dénomme 
habituellement une inversion de population. 

On peut également observer des températures négatives pour des substances paramagnétiques 
lorsque le couplage entre les spins des atomes et les autres degrés de liberté (vibrations des 
atomes du réseau, par exemple) est plus faible que celui entre deux spins voisins. Pour des 
substances magnétiques, dans lesquelles coexistent des régions avec des projections de spin 
préférentiellement dans l’une ou l’autre direction, on a, pour celles-ci, une température qui 
peut être positive ou négative. Ainsi, pour un même matériau, on peut trouver un mélange des 
régions de températures positives qui avoisinent des régions de températures négatives. 

Cette notion avec nous depuis 1950, Grâce à Ramsey Norman, et Malgré la réaction rapide 
de Norman Ramsey pour trouver dans les principes de la thermodynamique et ceux de la 
mécanique statistique une explication à l’existence des températures négatives, et montrer 
qu’elle n’entrait pas en contradiction avec ces principes mais demandait une simple 
reformulation, plusieurs chercheurs n’étaient pas satisfaits. Certains envisageaient même que 
l’on puisse peut-être faire revivre le concept de mouvement perpétuel avec les systèmes 
quantiques à température négative, malgré les arguments de Ramsey indiquant qu’il n’en était 
rien, ce qui poussait les chercheurs scientifiques à chercher dans cette merveilleuse notion. 

Et voilà un groupe Allemands de l’université Munich en collaboration avec l’institut Max 
Planck d’optique quantique viennent de couper court à ces spéculations dans une expérience 
à base d’un gaz atomes ultra froids, où ils ont arrivé à des températures négatives. 

Dans le cadre de notre projet de fin d’études, nous nous intéresserons à l’interprétation 
statistique de cette notion, en se basant sur quelques concepts utilisés dans la démonstration 
de son existence. 

Notre travail comporte 4 chapitres principaux ; un premier chapitre consacré à expliquer les 
deux sens physiques de températures, le sens thermodynamique où la température est décrite 
par l’agitation des particules constituant la matière, elle suit le théorème d’équipartition, alors 
que le sens statistique repose sur le lien entre entropie et énergie. 

Le deuxième chapitre abordera la notion d’entropie et ses interprétations, statistique où elle 
exprime le manque d’infonnation, et thermodynamique où elle est décrite par les 2 principes 
le second qui nous donne dans quel sens la transformation faite et le 3ème de Nernst où 
l’entropie à 0 kelvin est nulle. 
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Le troisième chapitre porté sur la notion de la température négative, son historique, sa 
définition, démonstration de son existence théoriquement. 

Le quatrième chapitre repose sur, le processus d’inversion de population dans les lasers et par 
un champ magnétique, les définitions des notions dont basée l’expérience, et une description 
de l’expérience. 

Le cinquième chapitre a pour but de créer un programme en langage c pour étudier la 
variation de l’entropie en fonction de l’énergie, à l’aide des deux ensembles: microcanonique 
et canonique, ainsi de voir l’influence de la répartition des niveaux sur les résultats. 

Le sixième chapitre est consacré pour le champ d’application de cette notion. 
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CHAPITRE 1 : Définition générale 


I. La température : 

La température mesure une grandeur physique liée à la sensation de chaud et froid d'un corps 
par rapport à un autre corps. Par définition, une température est une mesure numérique d'une 
chaleur, sa détermination se fait par détection de rayonnement thermique, la vitesse des 
particules, l'énergie cinétique, ou par le comportement de la masse d'un matériau 
thermométrique .La température d’un corps est la mesure de son agitation moléculaire, c’est 
une fonction de la proportion d’atomes excités et d’atomes à l’état fondamentale d’un corps 
.Elle représente l’énergie moyenne d’une molécule . L'unité légale de température dans 
le système international est le kelvin de symbole K (noter l'absence du symbole 0 car ce n'est 
pas une échelle de mesure). Il existe d'autres systèmes de mesures antérieurs et toujours 
utilisés : les échelles Celsius, centigrade, Fahrenheit et Rankine. [ - ] 

II. Ses types : 

Il existe deux types de la température la 1ère thermodynamique et la 2ème statistique : 

1. Interprétation thermodynamique de La température : [4] 

a. Définition : 

La température thermodynamique, usuellement notée T, 
est une mesure absolue de la température et constitue 
l’une des grandeurs principales de la 
thermodynamique. 

La température thermodynamique est une mesure 
absolue parce qu’elle traduit directement le phénomène 
physique fondamental qui sous-tend la température : 
l’agitation des particules constituant la matière. Son 
point origine, ou zéro absolu correspond à l’état de la 

matière où ces particules ont une agitation minimale et ne peuvent plus être refroidies. 
Simplement dit, la température traduit l’énergie cinétique des constituants de la matière 
(atomes et molécules). L’agitation cinétique, l’énergie potentielle de particules ainsi que 
certains autres types d’énergie de nature quantique emmagasinent l’énergie interne dans la 
matière : l’énergie interne totale que renferme un corps y est répartie entre diverses formes 
physiques de stockage, mais seule l’énergie cinétique se traduit sous forme de température. 

b. Théorème d’équipartition : 

•Le premier principe de la thermodynamique impose que l’énergie totale d’un système, 
somme de toutes les énergies cinétiques et toutes les énergies potentielles restent toujours 
constantes quoiqu’il arrive. 
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Ce principe reste parfaitement valable en mécanique quantique contrairement au principe 
d'équipartition de l'énergie (Valable seulement dans la limite des faibles densités et des hautes 
températures). 

•En thermodynamique, la température est définie comme une quantité proportionnelle à 
l’énergie cinétique moyenne d’une particule \E — A xT 


•Pour chaque degré de liberté la constante A est égale à ^ k , où k est la constante de 
Boltzmann, de valeur numérique 1,38 x 10 ~ 23 ]K~ 1 , 

3 3 

Donc à trois dimensions : E — -kT et s’il y a N particules dans le système : E — - NkT .S 

•l'énergie d'un système est emmagasinée sous forme d'énergie cinétique et d'énergie 
potentielle : 

-énergie cinétique de translation : < V2 mv 2 > — V2 kT, < V2 mVy > — V2 kT, ... 

-énergie cinétique de rotation :< V2 I œ 2 > — V2 kT, 

-énergie potentielle de vibration : < V2 K x 2 > — V2 kT, ... 

•Enoncé du théorème : 

« Si l’énergie d’une particule est une somme de tennes quadratiques des dégrés de liberté, la 
contribution de chaque terme à son énergie moyenne est - kT quelque soit la nature du terme». 


•exemples 


< 

E 

>- 

É . 

2 m 

+ 

K| 

2 

donc 

< 

E 

>- 

2-kT = kT 

2 

< 

E 

>- 

Pi 

2m 

+ 

K l 
2 

donc 

< 

E 

>= 

6-kT = 3kT 

2 


►Energie interne d’un gaz parfait : 

•Energie interne d’un gaz parfait (monoatomique) : 


N 

V 1 1 -^2 .1 

E = 2^ 2 mVi = N (-rnV 2 ) 

i=l 


3 3 

E = -NKT = -nRT 
2 2 


avec NK — nN a K — nR 


•Cas général : 


IN 

E = Ej = N (E) 

i=l 

avec (E) est l'énergie totale d'une molécule 
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• Cas particuliers. Nombre v de termes 


Type de molécule Expression énergie E quadratiques dans E 


monoatomique 

E = 'mv 2 +(mvj + ' mv" 

z * z y z z 

v = 3 

diatomique, rigide 

E = (mv; +*mv 2 + (mv z +U,ui( + U,o) 2 

10 

11 

?► 

autre, rigide 

E =(mv; +jmvj +(mv; -t- ^1,00“ +U 2 to; +(I,toj 

v =6 

diatomique, vibrante 

E = 'mv; + ( mv 2 +(mvj +U,uf +(L t0 2 +ik(r-r 0 ) 2 

v = 7 


m = masse totale de la molécule k v ^ nombre de 

( v x ,v v z ) = vitesse du centre de masse degrés de li5erté » 

I,, I 2 ,1 3 = moments d’inertie principaux 
(û)j, oi 2 , co 3 ) = vitesse angulaire de rotation dans le repère d’inertie 
cas diatomique : I 3 = 0, r = distance entre les deux atomes, u = masse réduite 


•Remarque : 

E 

E — AT -> T = - , la constante A est positive, la même chose pour l’énergie E (car E est une 
somme des termes quadratiques) Donc la température T est toujours positive. 


c. Conclusion : 

L’interprétation thermodynamique de la température ne donne aucune chance à l’apparition 
des températures négatives car T=E/A est toujours positive. 


2. Interprétation statistique de La température: Entrop , ie/ 

- c - 1 --— desordre 

En statistique, la définition la plus rigoureuse de la 
température repose sur le lien entre entropie et énergie. 

Lorsqu’on augmente l’énergie d’un système, on modifie 
également le nombre d’états dans lesquels il peut se 
trouver. La température est définie comme la variation de 
l’énergie avec le désordre. Si un peu d’énergie 
supplémentaire augmente beaucoup le désordre du 
système, sa température est basse. S’il faut rajouter 
beaucoup d’énergie au système pour changer significativement son désordre, sa température 
est élevée. 



Energie 


Elle est donnée par : 


1 _ ÔS 
T ~ dE 


Avec E l’énergie du système et S son entropie. 
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CHAPITRE 2 : L’entropie 



I- Historique : 

L’entropie est la dernière et la plus mystérieuse des cinq grandeurs physiques (température, 
pression, volume, énergie interne, entropie) définissant l’état d’un système thermodynamique, 
c'est-à-dire d’un ensemble matériel délimité capable d’échanger de la chaleur et du travail 
avec le milieu extérieur. 

L’entropie caractérise l’aptitude de l’énergie contenue dans un système à fournir du travail, ou 
plutôt son incapacité à le faire : plus cette grandeur est élevée, plus l’énergie est dispersée, 
homogénéisée et donc moins utilisable (pour produire des effets mécaniques organisés). 
Lorsqu’un système thermodynamique est isolé sans échange possible avec l’extérieur, il ne 
peut qu’évoluer spontanément vers le maximum de son entropie pour tendre vers un état 
d’équilibre définitif (2 e principe) alors que son énergie interne reste conservée (1 er principe). 

La notion d’entropie a été introduite par Rudolf Clausius en 1865. Il avait entrepris d’intégrer 
dans une nouvelle discipline, la thermodynamique, l’ensemble des progrès faits en physique 
de la chaleur. Clausius s’attaqua à la traduction mathématique de l’équivalence chaleur-travail 
de Joule (W=JQ). Grâce à la machine à vapeur, on savait quantifier le travail mécanique d’un 
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gaz poussant un piston dans un cylindre à pression constante (À W = p À V) mais cette 
relation n’avait pas encore été transposée à son équivalent, la chaleur. 

Clausius, assimilant la température à une pression, proposa alors par analogie la relation 
ÀQ=T AS, S étant une grandeur d’état additive (extensive) nouvelle qu’il baptisa entropie 
(« transformation » en grec), son accroissement s’exprimant en Joules par degrés Kelvin 
(J/K). A la lumière du second principe imposant à la chaleur d’aller toujours du chaud vers le 
froid, la variation d’entropie s’est alors révélée être la grandeur traduisant le sens de toute 
transformation irréversible d’un système thermodynamique par l’inéquation AS>0. Mais à ce 
stade l’entropie macroscopique restait une abstraction, son sens physique restant une énigme. 

Clausius établit que la chaleur est un phénomène mécanique dû au mouvement des particules 
(molécules) et la température une mesure de leur agitation. Simultanément Maxwell, 
appliquant à la physique des particules les nouvelles méthodes de l’analyse statistique, parvint 
à retrouver la plupart des résultats macroscopiques de la thermodynamique, en particuliers la 
loi des gaz parfaits. Mais à l’échelle microscopique, ils ne purent mettre en évidence 
d’irréversibilités pouvant expliquer le deuxième principe et ou se cachait l’entropie. 



Ce fut Boltzmann qui découvrit en 1873 que la fonction de distribution statistique des 
positions et des vitesses des molécules obéissait à une évolution irréversible : pour un même 
état macroscopique (volume V, énergie interne U), constitué déjà d’un nombre gigantesque de 
particules (>10 23 dans une seule mole de gaz parfait), il existe un nombre Q de distributions 
possibles des positions et des vitesses encore beaucoup plus « immensément grand ». Or la 
probabilité des « équirépartis » dans tous les états possibles croît mathématiquement avec leur 
nombre au détriment des autres. A l’échelle microscopique, ces états O équirépartis sont si 
dominants en nombre et en temps d’occupation qu’ils constituent l’état d’équilibre 
macroscopique le plus probable vers lequel, par agitation thermique, le système a tendance à 
évoluer spontanément. Le ressort probabiliste du second principe qui pousse le chaud vers le 
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froid était découvert. Boltzmann a alors lié l’entropie d’un système au nombre Q par la 
relation : S = K x log il , où K est la constante de Boltzmann (1,38110.10' 23 ) 

On retrouve que plus les particules d’un système sont, en positions et en vitesses, dispersées 
de façon homogène entre leurs Q états possibles, plus l’entropie de ce système est grande et sa 
capacité à fournir du travail limitée. 

En 1904, le troisième principe (Nernst) stipulait que l’entropie est nulle au zéro absolu 
puisqu’un ordre parfait y règne, la dotant ainsi d’une échelle de mesure absolue. 

II- Interprétation statistique de l’entropie : 

Boltzmann définit l'entropie statistique d'un état macroscopique par la formule : 

S = K.logÇï 

Où il le nombre d'états microscopiques correspondant à l'état macroscopique considéré. La 
constante de Boltzmann est introduite afin d'identifier les températures cinétiques et 
thermodynamique. 

Ainsi définie, l'entropie est extensive, car, si l'on considère un système constitué de deux 
sous-systèmes indépendants, on a : 

H = a 2 et S = K. logïl = K. log + K. logïl 2 = + S 2 

Généralement cette entropie exprime le manque d’information pour cela : 

Considérons un ensemble de M événement de probabilité respective Pi . Pm ,dont la 

somme est égal à 1, lorsque les événements ont même probabilités, notre perplexité est 
d’autant plus grande que ces événements sont plus nombreux .le manque d’information qui va 

caractériser le caractère aléatoire des probabilité Pi . Pm est une certaine fonction 

S (Pi . P M ) de ses probabilités que l’on appelle l’entropie , ou l’incertitude ou le 

désorde associé à la loi de probabilité considérée. 

La démarche originale est celle de Shannon, à un facteur multiplicatif, K positif, l’entropie est 
définie comme un gain moyen d’information en 1949 par : 

S(Pl . P M ) = 2jm=l Pm logPm (A) 

Un cas particulier important est celui où les M événements sont équiprobables .L’entropie 
statistique se réduit alors à : 

Dm = s(s. i j )=K.logM (B) 

A et B sont introduises en 1949 dans le contexte de la théorie de la communication, en 
considérant chacune des messages m susceptibles d’être transmis par un émetteur à un 
récepteur comme un événement de probabilités ( P m ), donc la quantité de mouvement gagnée 
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lors de la récéption du messages m doit être définie comme une situation après récéption du 
message, par : 


Im= K. log(Prn) 

Les messages dont Pm est faible sont ceux qui transportent plus d’information ce qui conduit 
à une faible contribution à l’information 

carP m ->0 —logP m ->± oo, mais — P m logP m Q. 

l’entropie statistique caractérisera, pour un macro-état préparé selon une loi de probabilité 
donnée ,notre incertitudes sur l’ensemble de toutes les expériences concevables, et avec 
k=constante de Boltzmann sera : 


M 


S = 


- k l 


PmlogPr 


m =1 


1. Mesure de l’information manquante : 

- Maximum : Pour N fixe, S est maximale est égale à S — K. logN lorsque tous les 
événements sont équiprobables. = P 2 — = P N Donc S — k. log N 

- Minimum : Dans le cas d’un événement certain c.à.d. 5(0 .,1,0000000) = 0 

c.à.d. pi — 1, aucun manque d’information. 

- Croissance : Dans le cas d’un événement équiprobables l’entropie augmente avec N. 
- Evénement impossible : dont la probabilité est nulle : 


S(P 1 ,P 2 . P N , 0,0,0,0) = S(P 1 ,P 2 . P N ) . 

Symétrie : S(p 1 . p t ....pj . P N ) = S(p ± . pj ....p t . P N ) 

Additive : on regroupe les événements de 1.N 

Événement A: p 1 + p 2 +... +p N — q A et n+1 .N 

Événement B: +. l N — q B . Alors : 


S(Pv P 2 . P N )= S(q A , q B ) + q A -S(- . —) + q B .S .—) 

\qi q A ) \ q x q B J 

-Sous additivité : on considère un événement composé de 2 événements : 

-Zeler€a = {l. m, .M} , le2ème€b={l .n,.N} 

-Pmn,> probabilité des événements n :p“ = Y^^Pmm Pu = ün+iPmn 

-On peut montrer que les entropies statistiques satisfont à l’inégalité de sous -additivité : 
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s (Pi . Pm) + S(p? . Pn)> ^(Pll - ■ Pmn . Pmn ) 

-AL’égalité est atteinte pour des événements a et b qui sont indépendants : p mn — PmPn 

Concavité : Deux lois de propriétés :p x ., p N , p\ . P' N associés aux mêmes 

événements, on construit une nouvelle loi de probabilité : 

p" x = ap x + (1 — cï)p' 1 P" N — ap N + (1 — a)p' N avec 0 < a < 1 
-L’entropie statistique est alors une fonction concave des probabilités. 

S(p'\ . P" N ) > aS (.Pi . Pn) + (1 - a)S(p\ .pV) 

Remarque : S ——k'Zf=iPi^ogp i 

Théorie de l’information : k=l ~^Log — log 2 

En thermodynamique: k=constante de Boltzmann=l ,38/é = 1,3810 -23 /. k~ 1 

2. entropie statistique d’un état quantique : 

S = —fcüf =1 Pj log p, et p=£”=i Pi \iii >< Ui | dans cette base (où p est diagonale), 

S = -KT r (plogp) , C’est l’entropie de Von Newman. 

Lemme :lnégalité de l’entropie statistique : 

Quelque soit 2 opérateurs densité p et p’d’un même espace de Hilbert : S(p) < -kT r (plogp) 
L’égalité est atteinte si : P=P' 

Propriétés : 

- Maximum : si les kets possibles appartiennent à un sous espace de Hilbert H w de H de 
dimension w, l’entropie statistique est maximale et égale à S = klogw,si les probabilités de 
tous les kets sont égales . 

- Minimum : l’entropie statistique est nulle pour un système dans un pur. 

- Invariance : l’entropie statistique est invariante si on applique une transformation unitaire à 

U + pU 

- Additivité : a , b, 2 systèmes statistiquement indépendante décrits parp a et p b , 

P = P a x p b dans H = H a x H b ^S(p) = S(p a ) + S(p 6 ) . 

- Corrélation : sous additivité ; soit un système composé par 2 sous systèmes a et b, 

p€H = H a x H b S(p) — 5(p a ) + S(p b ), l’égalité est atteinte si : p = p a x p b 

- Concavité : l’entropie statistique est une fonction concave sur l’ensemble des opérateurs 
densité d’un espace de Hilbert H quelque soit p 1 et p 2 et 0 < a < 1 
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s(ap x + (1 — a)p 2 > aS^) + (1 — a)S(p 2 ).l’égalité est atteinte lorsque p 2 = pi- 


Evolution temporelle de l’entropie : 


Pour un système libre : * [ZZ, p] , si H est parfaitement connu : 


ds 

ât 


= -kTr(^logp + ^j = -kTr^(log p + 1) = -klog^logp 


ds 

— — —kTrH\ p, Zoop] = 0 -> s — cte 
ot 


3. entropie statistique d’un état classique : 


En classique l’entropie s’exprime en fonction de la densité en phase par : 

S(p) = -k£ N /dT N P N logP N (A) 

Si N est connu -> (A) se réduit à 1 seule intégrale sur 6N variable de l’espace des phases. 
Contrairement à la mécanique quantique l’entropie est non négatif c.à.d. n’est pas bornée 
inférieurement .par suite, le fait d’obtenir un s(p) négative signale que la densité en phase p ne 
peut être considérée comme une limite classique de matrice densité et que l’approximation 
classique n’est pas donc valable. 

Brièvement, la théorie de l’information permet d’associer à tout état d’un système représenté 
par un opérateur densité, un nombre, le manque d’information ou entropie statistique qui 
mesure le degré de désordre venant du caractère aléatoire de l’état .l’entropie statistique varie 
entre 0 (cas pur) et celle donnée par Boltzmann. Elle croit dans les circonstances où le 
désordre augmente, suppression des corrélations entre sous -systèmes, réunion de plusieurs 
ensemble statistique en une seule évolution imparfaitement connue .l’entropie de sous 
systèmes statistiquement indépendants s’ajoute.la limite classique est donnée par : 




dT N pj\flogp N 


III. Interprétation thermodynamique de l’entropie : 

En thermodynamique l’entropie est décrite par deux principales : 

1. Second principe : 

En thermodynamique le premier principe permet de faire le bilan des transferts d’énergie d’un 
système subissant des transformations quelconque, mais ne fais aucune différence entre les 
diverses forme sous lesquelles l’énergie se trouve échangée. Cela fait appel à l’introduction 
d’une nouvelle grandeur l’entropie à l’aide du deuxième principe. 
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Ce principe introduit la fonction d'état entropie : S, usuellement assimilée à la notion de 
désordre qui ne peut que croître au cours d'une transformation réelle. 

Enoncé : 

« L’entropie augmente lors d’une transformation irréversible c.à.d. 
augmentation du désordre». 

2. Troisième principe : 

le second principe ne permet de définir la fonction entropie qu’à une constant près, qu’il 
s’agisse d’un corps pur ou d’un système plus complexe: la différence d’entropie entre deux 
états est toujours définie sans ambigüité, elle peut être évaluée en connaissant l’équation 
d’état du corps pur et en déterminant sa chaleur spécifique C v ouC p . 

L’indétermination de la fonction S se répercute sur d’autres fonctions d’état telles que: 
L’énergie libre F — U — T S L’enthalpie libre G — H — T S 

Ceci a incite les physiciens à chercher un état de référence d’un caractère universel pour lever 
cette indétermination, d’où l’énoncé de 3ème : 

«l’entropie S de tous les corps purs réels, chimiquement homogènes ,tend vers une 
valeur limite finie S 0 — 0 lorsqu’on s’approche absolu». 


3. Variation de l’entropie en fonction de la température : 

L'entropie d'un corps pur est une fonction croissante de la température, cette croissance est 
discontinue à chaque température correspondant à un changement d'état physique du corps 

pur. Le symbole de cette grandeur d'état est , La dérivée de l'entropie par rapport à la 

, ds c p 

température est : — = 


S (J/K) 


Liquide 



Gaz 

_ 



AS#C W" (VL) 

Vaporisation 

AS = L /T ^ 

(vatf <vatf 






Fusion 

AS = L 

/T 


(fus) (tus) 


b s#c P( »., ln (VL) 





T (K) 


w 


T(fus) 


T(vap) 


Variation de l’entropie molaire d’un corps pur en fonction de T 


Le diagramme proposé ci-dessus correspond au cas de figure où l'on observe la succession des 
phases solide —*■ liquide —> gaz au cours du chauffage d'un corps pur, effectué à pression 
constante. Cela signifie que la pression choisie est supérieure à celle du point triple figurant 
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CHAPITRE 3 : La température négative 


I. Historique : 


NORMAN RAMSEY 


Dans les années 1950, Purcell, Pound et 
Ramsey envisagent cependant l’existence 
des températures inférieures au zéro 
absolu. Leur proposition repose sur une 
définition plus fine de la température, 
élaborée au cours du XIXe siècle et qui 
fait appel au concept d’entropie. 

Lorsqu’on décrit une assemblée de 
particules (par exemple l’air de la pièce), 
on utilise des grandeurs d’ensemble : la 
température, la pression, le volume, etc. Ces grandeurs ne décrivent pas le système de manière 
unique : si une des particules avait une position légèrement différente, la température, la 
pression et le volume de l’ensemble resteraient inchangés. Il en résulte un certain désordre : le 
système peut être dans un grand nombre d’états différents, mais équivalents. 




En novembre 2012, l’équipe d’immanuel Bloch à l’université Munich en collaboration avec 
l’institut Planck d’optique quantique de Garching, présente dans la prestigieuse revue Science 
des résultats expérimentaux attestant de température négative. L’expérience est réalisée à 
l’aide d’un gaz d’atomes ultra froids de potassium dans un réseau optique. 


IL Défintion : 

Une température négative en thermodynamique est une température exprimée par un 
nombre négatif de kelvins. Le signe moins apparaît lorsque la définition de la température 
utilisée dans la thermodynamique moderne est appliquée à certains états quantiques très 
particuliers, dont l'entropie diminue au fur et à mesure qu'on leur ajoute de l'énergie. Ces 
valeurs négatives ne sont pas à considérer comme étant « plus froides » que les températures 
positives et elles ne signifient pas que la température est passée à un moment quelconque par 
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le zéro absolu. « Ce dernier restant impossible à atteindre ». Pour passer d'une température 
positive à une température négative, le système passe en fait par des températures infinies, 
d'abord positives puis négatives. 

Lorsque l'on se réfère familièrement à des températures exprimées par des nombres négatifs 
sur les échelles de Celsius et Fahrenheit, celles-ci représentent des valeurs plus froides que le 
zéro de ces mêmes échelles mais sont néanmoins plus chaudes que le zéro absolu. En 
revanche, un système dont la température est réellement négative dans l'absolu sur l'échelle de 
Kelvin est plus chaud qu'un autre système dont la température est, elle, positive, ce qui est 
paradoxal. Si un système à température négative entre en contact avec un système à 
température positive, la chaleur se propage depuis le système à température négative vers 
celui à température positive (voir III. 2). 

Le compromis entre énergie et entropie est expliqué comme suit : 

Les systèmes dont la température est positive voient leur entropie augmenter quand on leur 
ajoute de l'énergie. À l'inverse, les systèmes dont la température est négative voient leur 
entropie diminuer quand on leur ajoute de l'énergie. 

La plupart des systèmes que l'on connait ne peuvent pas acquérir de température négative 
puisque leur entropie augmente toujours lorsqu'on leur ajoute de l'énergie. 

III. Interprétation : 

l.Températures négatives dans les systèmes à deux niveaux : [1] [2] [8] [9]x 

L’exemple d’un système paramagnétique avec 
j=l/2 (système à deux niveaux) nous permet de 
discuter d’une extension possible de la notion de 
température. 

Chacune des N particules du système peut prendre 
deux valeurs possibles de l’énergie ±e avec s — 
gg b j où g b est le magnéton de Bohr et g le facteur 
de Landé. 

Soit N + le nombre de particules d’énergie +s et AL 
celles d’énergie —s . On a évidemment 
N = N+ + AL (1) 

L’énergie totale du système vaut : 

E = (iV + -AL)£ (2) 

On peut tirer les expressions de N + et AL des équations (1) et (2) 

*- = >+-!) 

Le nombre Q d’états accessibles pour un système d’énergie E donnée, et possédant N 
particules (micro canonique) est : 
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N! 


N! 


[f(N +E/e)]l [| (N -E/e)]! 


(3) 


Q(E,N ) = Ciï + = Ciï_ = N+] N _ } _ 

Il y a en effet exactement N! façons différentes de dénombrer les particules, où cependant les 
permutations des particules dans 1 ‘état +e (N + \) ne changent pas la situation (ainsi que les 
permutations des particules dans l’état -e (AL!)). 

C {/ + : Le nombre de tirages aléatoires de N+ particules parmi N (sans remise). 

N _ : Le nombre de tirages aléatoires de N. particules parmi N (sans remise). 


C N 


I)] 


+ 


L’équation (3) permet de calculer l’entropie du système. 

S(E,N) = klogQ(E,N) 

= k[NlnN —N— i(N +|)ln [)(« + f)]+ i ( N +f)" ï ( N “ f) 1 " [Î( N 

i (N - |)] 

= k [N ln N-i(N+f) ln [i(jV + f)]-i(N-f) ln [i(N-f)]-N+i(N+f) + 
ï ( N - !)] 

S(E, N) = k [n In N - i (N + f) ln g (iV + f)] - \ (n - f) ln [> - f)]l (4) 


où nous avons utilisé la formule de Stirling puisque N + ,N_» 1. (ln N! = N ln N — N) 
L’équation (4) donne pour la température 


ï= PK= ff| N 


= k 




1 

rl 

( E\] 

-ln 

- 

N + —) 

[ 2£ 

[2 

\ EJ 


2 £ + 


2 £ 


■ln 


r 1 / Eyi 

1 

-IN - 

+ —— 

12 V EJ. 

2£ 


2 £ 


î("+r) 


K N - E/e 

= 2F ln( NTlÂ ) 


) d ’ où r-ïhSgr (5) 


^ 1 K , f N - E/e 

Donc: - =—ln (-— _- 

T 2s v N+E/£ y K L v N+E/£ 

Tant que E < 0 ; c’est-à-dire, qu’il y a plus de particules dans le niveau bas -s , nous aurons T 
> 0 comme d’habitude. Si l’on choisit la température pour variable, alors pour T oo les 
particules se répartissent de manière égale entre les niveaux haut et bas 
(N + — AL => E — 0).De cette façon nous ne sommes donc pas amenés à une situation 
nouvelle. Cependant il est possible d’exciter les particules d’un tel système par des 
mécanismes spéciaux de telle sorte que le niveau haut +£ soit plus occupé que le niveau bas. 
Dans les lasers ce mécanisme se nomme pompage ; on peut par exemple le faire optiquement 
en utilisant le rayonnement lumineux. Une condition nécessaire, cependant, dans ce cas les 
particules restent suffisamment longtemps sur le niveau haut (niveau métastable) ; si non une 
occupation élevée N + » AL ne poura pas être atteinte, un état pour lequel N + » AL se 
nomme inversion de population, puisqueE = (N + —N_)£ > 0, à une inversion de population 
correspondra une température négarive d ‘après l’équation (5). 

Pour mieux le comprendre, traçons S /(NK) d’après l’équation (4) en fonction de E /{Ne). 
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Figure 1 


Pour E /(Ne) = -1 toutes les particules sont dans le niveau bas, La pente |3oc dS/dE est infini 


et T=0. 


Dans le cas d’une répartition égale entre le niveaux haut et bas (chaque niveau d’énergie est 
également occupé N + — AL), nous avons E /(Ne)=0, l’entropie tend vers une valeur 
maximale S — Smax = Nk ln2 et la pente (1=0. 

Or, comme l’entropie est une mesure du désordre, alors le système tend vers une 
désaimantation totale ; un tel état est caractérisé par le fait qu’il n’y a aucune direction 
privilégiée des moments magnétiques des dipôles. 

En ce point (E=0) la température passe soudainement de +co à -oo (le passage des 
températures négatives à des températures positives se fait à l’infin), le système passe en fait 
par des températures infinies. 

Dès qu’il y a légèrement plus de particules sur le niveau haut que bas, 

Enfin si toutes les particules se trouvent sur le niveau haut (pompage), on a E /(Ns) — 1 et de 
nouveau T=0. 

On peut aussi voir ces températures négatives à partir de l’équation (5). 

Traçons T en fonction de E. 


T{E) 


E 


-Ne 


+ Ne 
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Il est donc raisonnable de décrire par une température négative un état un état inversé dans un 
tel système. De tels états inversés n’ont pas seulement de l’importance pour les lasers ; des 
substances possédant un moment magnétique nucléaire peuvent également subir une inversion 
de population Pour ce faire, il suffit d’aligner autant de dipôles que possible par un champ 
magnétique élevé à des températures pas trop élevées. 

Si maintenant on inverse brusquement l’orientation du champ, presque toutes les particules 
vont se trouver dans la «mauvaise» orientation (inversion). Elles vont évidemment après 
un certain temps de relaxation se remettre suivant la «bonne» orientation, mais pendant ce 
temps in peut faire correspondre une température négative. 

Pratiquement cela s’observe pour le fluorure de lithium (Purcell et Pound 1951). Les spins 
nucléaires sont alignés par un champ magnétique intense. Un équilibre thermique des spins 
entre eux est atteint au bout d’un temps de relaxation n. Qui est de l’ordre de 10 5 s. Si 
maintenant on inverse brusquement le champ, on a un état inversé dans lequel les spins son 
effectivement en équilibre entre eux au bout de 10‘ 5 s. Il est alors raisonnable de faire 
correspondre une température négative aux spins dans cet état. Il faut faire attention au fait 
que tout cristal de LiF se trouve encore à la température du laboratoire, seul le système des 
spins est influencé par l’inversion du champ. Le temps de relaxation nécessaire à l’obtention 
de l’équilibre thermique entre les spins et le réseau cristallin est alors très grand, T2=5 min. 

Par conséquent, pendant le temps T 2 il est possible de déterminer la température négative du 
système de spin par le moment magnétique nucléaire négatif. Ce n’est que 5 minutes plus tard 
que le système se retrouve en équilibre à la température de laboratoire. 

Remarquons enfin que ce n’est que pour les systèmes possédant une borne supérieure de 
l’énergie que des températures négatives peuvent être raisonnablement définies. Ce n’est 
qu’alors que la courbe S(E) se comporte comme à la figure 1. 

Les systèmes de spins ont la particularité que leur spectre d'énergie est borné supérieurement 
(contrairement à tous les systèmes possédant des degrés de liberté d'énergie cinétique), ce qui 
permet de définir pour eux des états à la température négative. 

Un système à température négative est un système qui, contrairement aux systèmes à 
température positive, a d’autant plus de chances de se trouver dans un état d’énergie donné 
que cette énergie est plus élevée. On comprend que pour un système «normal» dont le 
spectre d’énergie n’est pas borné supérieurement, les températures négatives soient 
irréalisables car elles conduisent à une énergie moyenne infinie, mais elles sont non seulement 
concevables mais ont été réalisées et étudiées pour les systèmes de spins nucléaires. Par 
ailleurs il faut voir qu’un système à température négative est plus chaut que n’importe quel 
système à température positive et que, mis en contact avec un système «normal» dont le 
spectre d’énergie n’est pas borné et donc nécessairement à température positive, il lui cédera 
de l’énergie de façon irréversible jusqu’à ce que tous ses états soient également peuplés, ce 
qui correspond comme on sait à une température infinie. Ce n’est qu’ensuite qu’il pourra se 
mettre en équilibre thermique avec le système normal à une température positive. Le passage 
pour ce système par une température infinie c’est-à-dire un état de chaos total est un passage 
obligé. 
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2. Échange de chaleur entre deux corps : [2] 

le but de cette partie est de montrer qu’un système à température négative est plus chaud que 
n’importe quel système à température positive. 

On considère deux corps formant ensemble un système isolé, mais ne se trouvant pas en 
équilibre mutuel, leur température Ti et T 2 sont différentes, l’équilibre s’établira dans le temps 
entre les corps et peu à peu les températures s’égaliseront. 


En applique le principe de Carnot : — > 0 => — + — > 0 (6) avecS=Si+S 2 
On sait que : E, = TjSj + pjVj , i = 1,2 
Donc dEi — TidSi avec dvi — 0 (vi est cte ) 


1 rr, 1 dEi . 1 dEl 1 dE2 ^ _ 

cLEj — TidS; => dS; = — — donc d apres (6) : —-—— —— > 0 (7) 

1 11 1 Ti dt v v y Tl dt T2 dt v y 

or l’énergie totale est constante E — cte — E ± + E 2 donc dE — 0 = dE t + dE z 

D’où dE 2 — —dEi donc (7) devient Tz Tl ËËl > q 

z 1 v T 2-Ti dt 


T2-T1 n dEi n 

Si Ti<0 et T2>0 : ——— < U —— < U donc le système à température négative cède de 

l'énergie sous forme de chaleur au système à température positive. 


CHAUD 


FROID 


-► 

Échange de chaleur 

D’après, la formulation de Clausius du second principe qui veut que la chaleur ne passe 
spontanément que d’un corps chaud à un corps froid. Donc le système à température négative 
est plus chaud que le système à température positive. 


24 














CHAPITRE 4 : L’expérience 

De 1950 jusqu’à 2012, des chercheurs scientifiques pensent à l’expérience pour montrer 
l’existence de ce type de température, il y a ceux qui ont pris une année sabbatique pour cette 
raison, comme David Auerbach en 1992, mais il n’a rien trouvé. 

En 2013, un groupe Allemand de l’université Munich en collaboration avec l’institut Planck 
ont démontré cette existence à l’aide des atomes « ultra froids » et des lasers à vides poussés, 
en se basant sur la notion d’inversion de population (déjà parlé dans l’interprétation), spins 
nucléaires. 


Tout d’abord qui ce que ça veut dire les lasers à vides poussés, les atomes ultras froids ? 

I- Inversion de population : 


Pour la compréhension plus profonde de la 
température négative, nous devons tous d’abord 
saisir pleinement ce qu’est une inversion de 
population. 

Une inversion de population survient lorsque 
des particules sont dans un état d’énergie plus 
élevé qu’un même système dans des conditions 
similaires, par exemple un système (à trois 
niveaux d'énergie minimum) comme un groupe 
d'atomes ou molécules, se trouve dans un état 
dans lequel la majorité des éléments sont dans 
un état excité plutôt que dans leur état fondamental N 2 >Ni. 


Energie 


„ Inversion de population 
entre ces deux niveaux 


Niveau fondamental 


Population 


Avoir une population plus élevée dans le 
niveau de plus haute énergie E 2 que dans le 
niveau Ei n'est pas une situation d'équilibre. 
En effet, à l'équilibre thermodynamique, la 
répartition des populations sur les niveaux 
est donnée par la loi de Boltzmann : 


N 2 


— N t ex p 


E 2 ~E t \ 
K b T ) 



Dans ce cas, N 2 est toujours inférieur à Ni. Il faut donc créer une situation hors équilibre en 
apportant de l'énergie au système des atomes via un pompage dont l'objectif est d'amener 
suffisamment d'atomes sur le niveau du haut. C'est ce phénomène d'inversion de population 
qui conduit à des températures négatives (en kelvin). 


25 


exosup.com 


page facebook 














En produisant une inversion des populations d’atomes entre deux niveaux, comme on le fait 
dans les lasers, on réalise une situation instable que l’on peut caractériser par T < 0. 


a) population normale 
(équilibre thermique) 

M 3 -- 

N 2 -- 

M-, -" »■ ■ » > 

N 1 > NI 2 > N 3 

Le facteur de Boltzmann : 

Étant donné un système en équilibre à la température uniforme T, la probabilité Pi pour une 
particule d’avoir l’énergie Ei est proportionnelle au facteur de Boltzmann : Exp(—Ei/kBT) 

Pi oc exp(—Ei/k B T) 

Et d’après cette formule, Les niveaux d’'energie plus élevée sont donc les moins peuplés 
lorsque T > 0. 


b) population inversée 


:n ergy 
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■ 2 
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2 
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m 3 

m 2 
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II. Inversion de population dans un laser : 

1. Définition d’un laser : 

LASER est l’acronyme de Light Amplification by Stimulated Emission of Radiations qui 
peut se traduire par Amplification de lumière par émission stimulée de radiations. Le 
développement de l’acronyme fait apparaitre le mécanisme principal sur lequel repose 
l’émission laser, l’émission stimulée qui a été imaginé en 1917 par A. Einstein. Il faudra 
pourtant attendre 1960 pour qu’une émission laser soit démontrée. 

2. Principe d’émission stimulée : 

Le laser est émis par le principe d'émission stimulée. 
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Une onde laser peut être émise selon trois principes. Comme pour une onde lumineuse 
classique, on retrouve les phénomènes de l’absorbation et de l’émission spontanée. 

• Dans le cas de l’absorbation, un atome qui se trouve dans un état stable absorbe un 
photon et se trouve ainsi dans un état excité E 2 , sur un autre niveau. Dans ce cas aucune 
lumière n’est émise : le photon est absorbé. 

• Lors de l’émission spontanée l’atome absorbe de l’énergie et passe également d’un état 
stable E] à un état excité E 2 . Lorsqu’il retourne à son état stable, il libère de l’énergie sous 
forme de photon. 

• Avec le principe d’émission stimulée l’atome est reéxcité avant d’être désexcité : il 
libère alors deux photons qui auront la même fréquence, la même énergie, qui se 
propageront dans la même direction et vibreront en phase. 


3. Principe de fonctionnement du laser : 

Comment fonctionne un LASER ? 


source d’énergie extérieure 



Dans le milieu amplificateur (ou mécanisme de pompage) on trouve des atomes que l'on va 
exciter grâce a une source d'énergie extérieure (elle peut s'agir d'énergie électrique : décharges 
électriques dans un gaz ionisé, d'énergie chimique : une réaction dite "exo énergétique" libère 
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de l'énergie, ou enfin d'énergie lumineuse). On appelle procédé de pompage optique ce 
mécanisme, il fut mis au point par Kastler en 1950 ce qui lui valut le prix Nobel. 


L'émission stimulée « photocopie » les photons 



2 X h V 



Amplification par l’émission stimulée 


^X/X/Xj^ 


✓x/xA^r 
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■'X/X/X jy 
AAAy 
'A/X/X jy 
''xy\/A r 

-'xy\y\j r 

^x/x/v^ 


Ces atomes excités vont émettre des photons qui effectueront des allés et retours entre les 
deux miroirs, ces photons vont servir à l'émission stimulée d'autres photons et ainsi de suite 
pour créer un faisceau laser intense. 

Le miroir de sortie est un miroir réfléchissant à 99% (partiellement réfléchissant) percé en son 
centre de façon a ne laisser s'échapper qu'un mince faisceau lumineux laser. 

Quand la source d'énergie excite en continu les atomes la lumière laser produite est dite 
continue et aura une puissance constante. 



4. Processus d’inversion de population dans un laser : 

Le phénomène d’inversion de population est une étape nécessaire dans le fonctionnement 
d’un laser. Dans cet appareil, l’inversion de population est un peu plus compliquée. En effet, 
nous avons des atomes à un niveau d’énergie initiale minimale (appelons E0). Pour faire 
augmenter l’énergie d’une molécule, nous devons fournir un photon d’énergie égale à la 
différence entre le niveau excité (appelons El) souhaité et celui de départ. En temps normal, 
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l’électron ainsi excité redescend très rapidement (le temps de vie moyen sur un niveau excité 
est d’environs 10 7 s) à l’état initial en réémettant un photon identique. 

Cependant, il peut exister dans certains matériaux, un état (E’ état métastable) intennédiaire à 
Eo et Ei ou le temps de vie est beaucoup plus long (10" 3 s). Une durée de vie plus longue sur ce 
niveau aura pour conséquence de provoquer une accumulation d’électrons. 

Lorsque les électrons redescendent de cet état métastable vers l’état fondamental du système, 
il se produit l’émission de photons cohérents et en phases, phénomène baptisé LASER. 

Une schématisation de ce phénomène est présentée à la Figure 2. L’accumulation d’électrons 
sur cette couche plus élevée que EO ce qu’on appelle : L’inversion de population responsable 
de l’émission stimulée. 



groond 

State 

Figure2 : inversion de population dans le laser 


Puisque nous visons ici à donner des exemples d’inversion de population, nous attarder 
d’avantage sur le sujet du laser n’est pas justifié. 


III. Inversion de population par un champ magnétique : 

Nous avons, jusqu’à maintenant, donné des exemples d’inversions de population sans y voir 
aucun lien avec une éventuelle température négative. Pour faire cette transition, et parler 
d’une température négative, nous allons devoir faire l’étude d’un système physique à spins 
nucléaires dans lequel l’inversion de population est possible. 

1. Spins nucléaires : 
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charge 


spin 


Le spin est, en physique quantique, une des propriétés des particules, au même titre que 
la masse ou la charge électrique. Comme d'autres observables quantiques, sa mesure donne 
des valeurs discrètes et est soumise au principe d'incertitude. C'est la seule observable 
quantique qui ne présente pas d'équivalent classique, contrairement, par exemple, à 
la position, l'impulsion ou l'énergie d'une particule. 

Les particules sont classées selon la valeur de leur nombre quantique de spin les bosons de 
spin entier ou nul et les fermions de spin demi-entier. Fermions et bosons se comportent 
différemment dans des systèmes comprenant plusieurs particules identiques le comportement 
fermionque de l'électron est ainsi la cause du principe d'exclusion de Pauli et des irrégularités 
de la table périodique des éléments. L'interaction spin-orbite conduit à la structure fine 
du spectre atomique. 

Le système à spin nucléaire est utilisé pour la facilité de déterminer un micro -état, ce dernier 
est déterminé par la donnée d’un ensemble complet de nombres quantiques qui permettent de 
préciser de manière unique la configuration du système. 


2. Interprétation : 

Considérons un champ magnétique intense et constant dans l’espace et le temps. Considérons 
aussi une tige de matériel paramagnétique positionnée dans ce champ magnétique. 

De par sa nature, les dipôles du matériel s’orienteront tous de manière à se retrouver parallèle 
au champ magnétique. La Figurel représente cet état. 

Dans cette orientation, le système est à l’équilibre en étant à son état énergétique le plus 


faible. 


Si nous inversons spontanément le champ magnétique de manière à retrouver le nord à droite 
et le sud à gauche, nous aurons, pendant une certaine période de temps, les dipôles 
moléculaires orientés antiparallèlement au champ. Dans cette configuration, le système est 
dans un état d’énergie plus élevé que celui dans la Figure al et donc, nous avons une 
inversion de population (Figure a2). 

Notez qu’habituellement, le temps requis pour réorienter les dipôles est très court (en fait, il 
est dans la plus part des cas non mesurable expérimentalement).Et donc le système revient 
rapidement à son état normal (température positive). 
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Les dipôles ne peuvent se réorienter instantanément après l’inversion du champ magnétique 
externe. Donc, dans l’intervalle de temps entre le moment ou nous avons inversé le champ 
magnétique et le moment où les dipôles se réorganisent la température devient négative . 

Le temps de réorganisation des dipôles, qu’on appelle temps de relaxation, est généralement 
très court (de l’ordre de 10-5 secondes pour la plupart des matériaux). 

Cependant, en 1951, une équipe de chercheurs (Purcell et Pound) [ -' ont réussi à obtenir une 
inversion de population des spins nucléaires dans des cristaux ultra pure de LIF de l’ordre de 
5 minutes. Cet état quasi stable a permit de mettre en évidence la réalité des inversions de 
populations ainsi que des températures négatives. 



Figure al Figure a2 

État initial du matériel paramagnétique Système tout de suite 

dans un champ externe. Après inversion du champ. 

Avant de revenir à un état stable (parallèle au champ), les dipôles passeront par un état de 
désordre complet. La Figure a3 nous illustre ce passage et la Figure a4 nous montre le 
système à son état d’équilibre final, le même qu’au départ avec les pôles inversés. 



Figure a3 : 

Passage par un désordre complet 


Figure a4 : 

État d’équilibre final. 


Une température négative ne s’obtient pas en refroidissant une substance sous les 0 K, mais 
plutôt en l’amenant au delà d’une température absolue infinie (Figure a5). 
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Figure a5 : Interprétation conceptuel à la température négative. 


IV-L’expérience : 

1- Des physiciens allemands ont porté un gaz d'atomes de potassium à une température 
absolue négative. Qu'est-ce que cela implique ? 

En thermodynamique, on décrit la température d'un gaz d'atomes par l'énergie cinétique - ou 
agitation - moyenne des particules qui le composent. Ulrich Schneider et ses collègues de 
l'institut Max Planck pour l'optique quantique, à Garching, et de l'université de Munich ont 
démarré leur expérience avec un gaz d'atomes de potassium à quelques milliardièmes de 
kelvins au-dessus du zéro absolu. La plupart des particules qui le composent se trouvaient 
ainsi dans des états d'énergie très faibles, minimaux même ; quelques rares particules 
occupaient néanmoins un état d'énergie plus élevé. Ce qu'ont fait les physiciens allemands, 
c'est qu'ils ont inversé cette distribution. Ainsi, ils ont placé plus d'atomes dans les états 
d'énergie élevés que dans les états d'énergie faibles. Le système, qui possédait au départ une 
température positive, a alors atteint ce qu'on appelle une température absolue négative 

2- Comment sont-ils parvenus à leur résultat ? 

D'ordinaire, l'énergie cinétique d'un atome n'a pas de valeur maximale imposée : on peut en 
principe l'augmenter indéfiniment. Pour donner un sens à une température absolue négative, 
Ulrich Schneider et ses collègues avaient donc besoin de modifier la nature de l'énergie 
cinétique, en lui imposant un plafond ; avec un laser, ils ont créé à cet effet un réseau optique. 
Puis, en appliquant brusquement un champ magnétique, ils ont réussi à faire passer les atomes 
de potassium des états de basse énergie à des états de haute énergie. Le système ainsi 
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CHAPITRE 5 : simulation numérique 


I. Introduction sur les méthodes numérique utilisée : 

Dans cette partie on a fait appel à l’analyse numérique pour les calculs des dérives en 
utilisant la méthode des différences finies, a pour principe: 

Le principe des MDF consiste à approcher une solution d’une équation aux dérives ordinaires 
ou partielles en un nombre fini de points seulement. Pour cela, on approche alors les dérives 
successive par des accroissements de la fonction entre 2 points, appelés aussi différences 
finies. 

Cette méthode donne 3 formules pour le calcul d’une dérivés, où h est le maillage: 
a- la formule progressive :/ (x) = ^ x+h ^ _ f^ 2 \x) x ^ 

b- la formule régressive : / (x) = ^^ + /®(x) x ^ 

c- la formule centrée : /(x) = f( x+h ^^ x — — f( 3 \x) x 

cette dernière est beaucoup plus précise que la régressive et la progressive , et c’est elle que 
l’on a utilisée dans notre simulation numérique. 


/ 

II. Etude d'un système à l'aide de l'ensemble microcanonique : 


1. Système à p niveaux : 

Imaginant des particules, chacune avec un schéma des états suivants : 

On analyse N particules, ou Ni sont dans l’état £i et N 2 sont dans l’état 

£2,... 

N=Ni+N 2 +...+N p 


e 


p 



E=Nl£l+N 2 £2+.. .+Np£p 

Au total le nombre de complexions possibles est égale à : 

no? mv_ r N r N ~ N ^ r N ~ N ^~ N ^ r *i-Ni-N 2 - Np _ ± 

IZfh, N) - L Ni .L N2 . L Ns ...L Np 

_ r N r N-N 2 r N-N 1 -N 2 r N-Ni_-N 2 -W p _! 

~ L N 2 - L N ± ' L N 3 "■ L N„ 


N! (N-NJI (W-jV 1 -N 2 )! 


N^iN-Nj] N 2 KN-N 1 -N 2 y. A^lQV-tfi-Nz-Na)! 


(W-W 1 -N 2 -Np- 0 ! 

NplQv-Ni-Nz-Na-Np-i-Np)! 


N! 

Q(E,N) = avec i — 1, ...,p 

lli ™O 


Cette équation permet de calculer l’entropie du système. 
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S(E, N) = KlogH(E,N) 


S(E, N) = K 


N.logN -^Nj.logNj - Ni 


( 8 ) 


= k[N ln N - N - N x ln N 1 +N 1 - N 2 ln N 2 + N 2 - N p _! In N p _ x + N p ^ - N p ln N p + N p ] 

Donc : 

S(E, N) = k[N ln N - N 1 ln N x - N 2 ln N 2 - N 3 ln N 3 -N p _! ln Np_ x - N p ln N p ] 


où nous avons utilisé la formule de Stirling puisque N v ..., N p » 1. (log N! = N log N — N). 

Puisque on ajuste 2 équations (N= Ni+N 2 +...+N p et E=Ni£i+N 2 £ 2 +...+N p £ p ) la contribution 
des particules ne peut être déterminé de manière exacte car on aura autant d'inconnues que 
d’équations (sauf dans le cas du système à 2 niveaux car on a 2 inconnues :N+ ,N — ) 


L’équation (8) donne pour la température 

^ = PK = ^ | N On sait pas exprimer S en fonction de E car on n’a pas trouvé une 

T u E 

relation entre les Nj et M (par contre), donc on ne peut pas calculer la température. 

On note que dans le cas d’un système à deux niveaux on n’a pas rencontrés ce problème 
(Voir l’équation (4) dans la partie III. a chapitre 2). 

2. Application Système à 3 niveaux (£ , 2£, 3s) 

Le nombre total de particules N: 

P 3 

K = 'Y d N t ='Y d N i = N 1 +N 2 + N 3 

i=1 i=1 

-> N = Ni + N 2 + N 3 (9) 

Avec N t le nombre de particules d’énergie£ , N 2 celles d’énergie 2£ et N 3 celles d’énergie 3£. 
N + , N 0 et N_ sont inconnus. 

L’énergie totale du système E: 

E 

E — N 1 .s + N 2 . 2 s + N 3 . 3£ = (iVi + 2 N 2 + 3 N 3 )s -» - = N 1 + 2N 2 + 3 N 3 

s 

M = Ni + 2N 2 + 3N 3 (10) avec M = | 

On a 2 équations [(9) et(10)] à 3 inconnues, il faut fixer l’un des inconnues pour calculer les 
autres, par exemple on fixe N j. 

3x(9) - (10) -> N 2 = 3N - M - 2N t (11) 

2x(9) - (10) -> N 3 = M - 2N + Ni (12) 
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Donc la contribution des particules est entièrement déterminée (selon la valeur de Ni) par la 
donnée de : M est N. 


Le nombre £2 d’états accessibles pour un système d’énergie E donnée, et possédant N 
particules (micro canonique) est : 

Q(E,N) =^n i (E,N) 

Avec Qj(E, N) est le nombre d’états accessibles de la contribution i. 

Donc il faut déterminer les contributions possibles pour chaque valeur de (E et M sont 
fixés) et calculer les £2j(E, N) de chaque contribution, par la suite faire la somme des £2j(E, N) 
Pour trouver le nombre totale d’états acccssibles£2(E, N) , et en utilisant la relation de 
Boltzmann Gibbs on aboutit à l’expression de l’entropie du système : 

S(E,N) = klogQ(E,N) = klogEiOj^N) (13) 

3. Exemple déterminer l’entropie d’un système à 3 niveaux avec N=3 et E=5e 
4 étapes nécessaires pour résoudre ce problème. 

-» 1- il faut déterminer les contributions possibles : 

Pour calculer Q il faut fixer /Vj. il faut tout d’abord déterminer les valeurs permises de N ± . 
Les valeurs possible qu’il peut prendre sont celles qui vérifient les contraintes : (9) et (10) 
On va tester tous les possibilités c.à.d. — 0,1, 2,..., N 
On remplace N et M par leurs valeurs dans (11) et (12) 




(U)-» 

N 2 — 4 — 2N ± (13) 



(12) -> 

N 3 = -l + N 1 (14) 

On commence le test : 


o 

II 

sT 

(13) -> 

N 2 - 4 (14) -» 

N 3 — — 1 or N 3 est un nombre positive 

donc N t 

= 0 n'est pas possible 


tH 

II 

sT 

(13) -* 

N z = 2 (14) -* 

o 

II 

«T 

Ai = 2 

(13) -> 

N 2 — 0 (14) -> 

tH 

II 

N 1 = 3 

(13) -> 

N 2 — —2 or N 2 

Est un nombre positive donc — 3 n'est pas 


possible. 

Donc les valeurs possibles de sont 1 et 2. 
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Et les contributions (A^; N 2 ',N 3 ) possibles sont : (1 ; 2 ; 0) et (2 ; 0 ; 1). 
On note : contribution 1 : (1 ; 2 ; 0) et contribution 2 : (2 ; 0 ; 1). 

-» 2- on calcule les £2; (E, N) : 

Les nombres £2i d’états accessibles de chaque contribution sont : 


Qi(E,N) 

£2 2 (e,n) 


N! 

3! 

N^.Nz! N 3 ! 

“ 1! 2! 0! 

N! 

3! 

N-lINzI N 3 \ 

_ 2!0! 1! 


-* 3- on calcule £2(E, N) : 

Le nombre £2 d’états accessibles est : 

2 

Q(E, N) = ^ Qi(E, N) = Q^E, N)+£2 2 (E, N) = 6 

i=l 


£2(E, N) = 6 


-* 4- on calcule S (E, N) : 

Enfin on calcule l’entropie du système en utilisant la relation de Boltzmann : 

S(E, N) = klog£2(E,N) = klogô 


Ces calculs peuvent être faits en utilisant le programme créer dans la partie IV- Progmmation 
en langage C : 
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le résultat enregistrer dans le fichier txt 

l'energie du système est: E=5e 
la contribution possible 1 : 

NI =1 N2=2 N3=0 

NI !=1 N2!=2 N3!=1 N1!N2!N3=2 
omega_1=3 

la contribution possible 2 : 

NI =2 N2=0 N3=1 

NI !=2 N2!=1 N3!=1 N1!N2!N3=2 
omega_2=3 

==> le nombre des contributions possibles est: 2 
==> le nombre oméga d'etats accessibe est: 6 

le résultat afficher sur la console 



4- Difficultés et Conclusion : 

Les difficultés qu’on a trouvé dans le cas d’un système à 3 niveaux décrit par l’ensemble 
microcanonique: 

38 








1- on ne peut pas exprimer Ni, N 2 et N 3 en fonction de M et s (comme dans le cas du système 
à 2 niveaux). 

2- pour chaque énergie E donnée, on a plusieurs contribution (par contre pour le cas 2 niveaux 
,pour chaque E. N+ et N. sont bien déterminés) 

3- la l ere difficulté a pour conséquence un problème au niveau du calcul de la température T ; 
on ne peut pas faire la dérivée de S par rapport à E ( ^ = -§ | N l’entropie S ne contient pas 
l’énergie E ). 


Conclusion : A l’aide de l’ensemble micro canonique, les calculs devient très longs pour 
décrire les systèmes dont les niveaux supérieur à 2, dans ce cas on va utiliser l’ensemble 
canonique . 


/ 

III- Etude d'un système à l'aide de l'ensemble canonique : 

1- Étude générale : 

L’énergie de la particule est : E m — £ m avec m un entier (indice du niveau) peut prendre des 
valeurs comprise entre p et q 

p <m < q 

p et q sont des entiers. 

Avec : 

p : indice du premier niveau, 
q : indice de derniers niveaux. 

-Les valeurs de m vont être déterminées par le programme après la saisir de n 

Dans le cas où les niveaux d’énergie sont symétriques par rapport à 0 : 

n — 1 n — 1 

p =-— etq = + — 

Ou on a fait l’approximation suivante : pour n très grand, on considère que les n paires sont 

. . , v . , n-1 

impaires (on ajoutant 1 a n), pour que le terme ±—^~ soit entier. 
n-1 : on a retranché le niveau fondamental (qui se trouve au milieu) 

Yi— 1 

—j— : On a divisé par 2 car les niveaux sont répartis de manière symétrique par rapport à l’état 
fondamental. 
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Dans le cas où les niveaux d’énergie sont positives : 
p — 0 et q — n—1 


Les grandeurs physiques : 

-la fonction de partition d’une seule particule s’écrit : 


r « = I 


z„ = 2_ i e 

m=p 


-la fonction de partition de N particules est : Z — Z a 

<? 

Z — ^ e~ N P £ n 

m=p 

-l’énergie totale : 


N 


dLogZ _ dLog Z a N dZ a N d e ~ pEr] 
h —-——— — —N - 


dp dp 

<? 

N 


Z a dP 


dp 


N Y?m= p de 

Z a dP 


N Y 

= zT2, 


P £ n 


m=p 


E 1 ^ 

~N~~ZZ 


P £ rt 


m=p 


NB :Les formules en rouge sont ceux qui ont utilisées dans le programme. 


dE 


-la chaleur spécifique à volume constant : C v = — = — — — = 




1 dE 
KT 2 dp 


N d 
Cv ~ ~KTÏ~dp 




-1 


[Z e "' ,E ’ n 

V £ m e-P £ ™ 


\m / 

m 


N 

~KT l 


Ui 


P £ r\ 


^£ m e P £m + \Ye- 


P £ n 


-y 

d/sLi 


P £ n 
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N 

~KT l 


-2 


-1 


P Cm J ^ e m e P E ™-\Ye P £ m\ ^ 


<? 2 p (3£ n 

c m c 


N 

C V - 2 


-2 


-1 


Z 6- "'" Z 




P^n 


Z 6- "'" Z 




2 p /?£n 


-l’entropie : 


dKTLogZ 
S =- —?— = K 


dT 


= K 


dLoq Z 

T —-1- Log Z 

-1 

T’^p(-E') + LogZ 


= K 


dLogZ dp 
T —dT + L °° Z 


S = K[pE + Log Z] -> - = N 

K 


pE 

— + Log Z a 


S BE 

— ~TT + Log Z a 


NK N 

La valeur maximale de l’entropie: 

S maximale <-> — = 0 -> KN — 1— + Log Z a ]= KN — — 0 -> p = 0 

dE dEIN a a \ N r 

Donc S = S max lorsque P — 0 quelque soit la répartition des niveaux (équidistants et 
non équidistants) 

Smax = lim K [pE + N Log Z a ] = KN lim log Z a = KN log n 

/> —>0 p —>0 


KN 


— logn 


Avec Z a — n lorsque P — 0 ; démonstration dans la partie suivante 


Les points particulièrs : (P — 0, +oo) 

a) P = 0 

-la fonction de partition 

<? 

P — 0 -> Z a — lim ^ e~ /}£m — lim fe~^ £ p + e -/Î£ (p +1 ) + —I- e~^M = 1 + 1 + — 1-1 
/?->o Z_i /?->o L } 

m=p 

= n 

Z a fl 
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-l’énergie totale : 

Q 

E v i \ 

P — 0 -* — — lim Z a _1 2^ s m e ~^ Em — ~ lim{£ p e~^ £ p + £( p+1 )e“^ Ê (p+i) + —I- Cqe - ^} 


m=p 


-~i £ p + £ (p+ 1) + - " + ^) 


n 


nE ^ 
~N ~ 2-i 


m—p 


Dans le cas des niveaux symétriques cette somme devient : 


nE 

~N 


n-1 

2 


— / Si — £ n-1 + S n—3 + • ■ • + 0 + ••• + £n-3 + En-1 

22 22 


n-1 


T - + E ^) + ( E -^ + c ^) 


3 + En- 3 ) + ”• + 0 — 0 + 0 + •■• + 0 


Car les niveaux sont symétriques donc a_ m — —£ m 


Et par conséquent : 


nE 

— = 0 

N 


-l’entropie : 


i ffSE \ 

B — 0 -> -= lim-1 - Loq Z n ] — limlogZ„ = logn 

p NK p^o\ N y a ) p^o & 5 

s = logn 


b) P — —oo -> b' — —oo 

-la fonction de partition 


P — 0 -> Z a — lim V e P £m — lim (e ^ £ p + e ^ £ (p+i) + —h e 

/?—>—oo / j R —»—co L J 


m-p 
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e P E[ i est le terme le plus petit 
Z n — lim e ^ £< ? 


-l’énergie totale : 


Cl 


B — —oo — = lim \ — ^ e^e P £m i = lim (—„— 

P ki „ . - 1 z m ( p ^_ œ \ -fie q ci j 


N /?->-oo I 


m-p 


n = £ « 


-l’entropie : 


p = ~ œ ^m = 4-» (t + 109 z “) = +,0 - 9 = ~ *+) 


s — 0 


P = +00 b — +00 


-la fonction de partition 


B — 0 -> Z a — lim V e P £m — lim (e P £ p + e ^ ê (p+i) + — h e 

fi —>-f- oo fi —>-i- oo L J 

m=p 


Z a - lim e Pe p 

fi —>-j- oo 


-l’énergie totale 


B — +oo -> — = lim 
N fi^+c 


\hl 




£r,e 


£ m e I e „fi £p -p 


/?£p 


m-p 


E 

N p 


-l’entropie 


B — +oo ->-— lim (— —I- LoqZ a ^\— lim [Be v + loqe P £ p) = lim {/?£„—/?£„) 

H NK fi^+oo V Al J fi^+ oo VP P « ) fi^_ oo^P P p P J 


s = 0 


2- Application : 

A- 1 er cas des niveaux non équidistants : 
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a. les niveaux d’énergie positive : E m = £ m > 0 
On considère le cas où E m — £ m — m 2 e avec e est une valeur donnée 

-Dans ce cas la fonction de partition d’une seule particule s’écrit : Z a — Y?m=o e _/?m2£ 


71-1 


*.- 2 > 

m =0 


-bm z 


Avec : b = (le 

-la fonction de partition de N particules est : Z — Z a 

-l’énergie totale : 


N 


E 1 
N 


n-l 


n-1 


— — m 2 e e P m2e = — m 2 e P m2£ 

Z„Z-J Zn Z—l 


m =0 


M 


n-l 

=-T 

z n Z-j 


m =0 


m 2 e bm2 


m=0 


Avec : M — — 

Ne 


-l’entropie 


S RE 

- —-h Loq Z a 

NK N y a 

s — bM -I- Log Z a 


Avec : 


s = 


NK 

Les points particulièrs : 

-la fonction de partition 

-l’énergie totale : 


/? = 0 -> = 0 


nE h= 


Z a Tl 


n-l n-l n-l 


n-l 


N 


£ =Z £,= Z £i= Z i2£=£ I/ 


1=0 1=0 1=0 


1=0 
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n-l 1 

i 2 — -n(n — 1)(2 n 

i =0 

n-1 

nM b=0 = ^ i 2 

i =0 

- 1) en utilisant le site: wolframalpha. com 

1 


M b=0 =-(n-lX2n-l) 

M b=0 :c’est la valeur de l’énergie réduite M lorsque b=0. 

-l’entropie : 



s = logn 

/3 = — 00 —> b — —oo 

-la fonction de partition 

Z a — lim e - ^ n-1 ) 2£ 

fi —>— 00 

-l’énergie totale : 

£ 

— — (n — l) 2 s 

N K J 

M — (n — l) 2 

-l’entropie : 

s — 0 



/? — +00 b — +00 

-la fonction de partition 

Z a — lim e~P° 2e 

+ oo 

Za = 1 

-l’énergie totale : 
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M = 0 


-l’entropie : 


5 — 0 


Remarque : 

1 

On se propose de déterminer la position de M b=0 — - (n — l)(2n — 1) par rapport au milieu 

6 

d’intervalle d’énergie. 


la valeur de l’énergie au milieu est : 


max+min _ 0 2 +(n-l) 2 _ (n-1) 2 


Pour cela on calcule œ — M b=0 


(n-l) 2 


Et selon le signe de œ on peut déduire la position de M b=0 . 

(O = - (n - l)(2n - 1) -= - (n - l)[(2n- 1) - 3(n- 1)] 

6 2 6 

1 ,1 

£o — — (n — 1) [2n — 1 — 3n + 3] = — (n — 1) [—n -1- 2] 

6 6 


m — — — (n — 1) (n — 2) 

6 

(72—1)2 (72—1)2 

Donc m < 0 <=> M b=0 - < 0 => M b=0 < 

n> 2 , donc œ — 0 pour n — 2 et œ < 0 pour n > 2 
Ou bien 


.. (n - l ) 2 

M b= o =--- pour n — 2 

(n — l ) 2 

M b=0 < --- pour n > 2 

Par conséquent la positions des points b — 0 et s — s max — log n vont être : 

-décalés vers la gauche pour n > 2 . 

-au milieu de l’intervalle de l’énergie. 


Voilà quelques résultats qu'on a trouver en utilisant le programme créer dans la partie IV- 
Progmmation en langage C : 

Exemple 1 : n=2 
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<—b=beta.e—> <—M=E/Ne—> 

-27.500 1.000000 0 

-28.000 0.99999999999930855000 

-27.995 0.99999999999930511000 

-27.990 0.99999999999930156000 

-27.985 0.99999999999929812000 

-27.980 0.99999999999929456000 

-27.975 0.99999999999929101000 

-27.970 0.99999999999928746000 

-27.965 0.99999999999928391000 

-27.960 0.99999999999928035000 

-27.955 0.99999999999927669000 

-27.950 0.99999999999927314000 

-27.945 0.99999999999926947000 

-27.940 0.99999999999926581000 

-27.935 0.99999999999926215000 

-27.930 0.99999999999925837000 

-27.925 0.99999999999925471000 

-27.920 0.99999999999925093000 


<—s=S/NK—> 

0.0000000000200550687168288280000000000000 

0.0000000000201509919861564410000000000000 

0.0000000000202504679691628550000000000000 

0.0000000000203463912384904690000000000000 

0.0000000000204458672214968830000000000000 

0.0000000000205453432045032970000000000000 

0.0000000000206448191875097110000000000000 

0.0000000000207442951705161250000000000000 

0.0000000000208437711535225390000000000000 

0.0000000000209432471365289530000000000000 

0.0000000000210427231195353670000000000000 

0.0000000000211457518162205820000000000000 

0.0000000000212487805129057960000000000000 

0.0000000000213518092095910110000000000000 

0.0000000000214548379062762250000000000000 

0.0000000000215578666029614400000000000000 

0.0000000000216644480133254550000000000000 


En utilisant ces valeurs, on trace la courbe qui représente s en fonction de M en à l’aide d’un 
logiciel graphique (Gnuplot dans ce cas). 


On obtient la courbe suivante : 



On remarque bien que la courbe est symétrique, car la position de s max (oit bien b — 0) se 
trouve au milieu(m = 0 ). 


1 


1 


M b=0 = ^ (n - l)( 2 n - 1) = - x 1 x 3 
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Vérification: 





M b=0 - 2 

(n — l) 2 1 1 

- 2 = 2 _ 2 ^ 

m = 0 => courbe symétrique 


Exemple 2 : n=5 

On refait la même chose pour n=5. 



on remarque bien que la position de s max (ou bien b — 0) est décalée vers la gauche. 
Vérification: 

M b=0 ~ ^( n ~ l)(2n — I)=gx4x9 

M b=0 = 6 
(n — l) 2 

<u = M b=0 ---= 6 — 8 

<n — — 2 <0 => décalage à gauche 

b. les niveaux d’énergie symétriques : positives pour m > 

0 et négatives pour m < 0 

On considère le cas où 
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£m | 


—m 2 s, m < 0 

m 2 e, m > 0 


-Dans ce cas la fonction de partition d’une seule particule s’écrit : 


n-1 
2 


n-1 

2 


n-1 

-1 2 


Z a = y e~ pEm = y e~P Em + y e~ pEm = y + y e~P m2£ - 

m=-^i ™=° m=-5=i ™=° 


n-1 

-1 2 

4 = X e " m2 + X e 

n-1 m =0 




n-1 n-1 

“2“ YT - 


n-1 

“2“ 


Z a = ^ e fcm2 + ^ e- fcm2 = 1 + ^ {e fcm2 + e- fcm2 } 


m= 1 m=0 

Avec \ b = fis 

-la fonction de partition de N particules est : Z = Z a 

-l’énergie totale : 


m= 1 


N 


£ 1 

iv 


n-1 

2 


= — Y = — Y . 

z « 4 . 


n-1 

2 


c p P £ m 

c m c 


n-1 
2 

y —m 2 e bm2 + m 2 e~ bm2 


m =1 


m=0 


n-1 

2 


n-1 

2 


M 


— — < ^ —m 2 e bm2 + ^ m 2 e bm2 


m =1 


m =0 


n-1 


M 


= — y { m 2 (e~ bm2 - e fcm2 )} 
Za Z—i 

m=1 


Avec : M — — 

Ne 


-l’entropie : 


5 RE 

-—- hLogZ a 

NK N y a 

s — bM + Log Z a 
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Avec : 


s = 


NK 


Les points particuliers : 


-la fonction de partition 


-l’énergie totale : 


/? = 0 -> = 0 


Z a Tl 


n—1 

2 


n—1 

~2r 


n—1 n—1 

“2r “2“ 


M - lim 

[3—> + CO Z a 


__ J ^ _ m 2 gbm 2 _|_ ^ 771.2 g bm 2 l — _ J _ ^ 777 2 + ^ 771 2 


m= 1 


m=0 


m=l m=0 


-l’entropie 


n-l n-1 

2 2 

7iM — — ^ 771 2 + ^ 771 2 
m=l m=l 

Mfj=o — 0 


s — logn 


P — — 00 —» è = — 00 


-la fonction de partition 


n-l 

-1 2 


n-l 

2 


Z a = lim 

f }->—OO 


y e èm2 + y e ~ bm2 = lim y 

/ I / I /?—> —OO / I 




n-l m=0 


m=0 


Za = lim é 

[3—>—00 




-l’énergie totale 
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M = lim 

/?->+oo 



-1 



g bm 2 


n—1 

2 



771 = 0 


,-bm 2 


— lim 

£>-> + oo 



,bm z 


M = lim 

b->+ oo 




M — — 



2 


-l’entropie : 


.v = 0 


Remarque : 

La courbe de s=f(M) est symétrique par rapport à 0, car eu = 0 : 


O) = M b= o - 


M-min "1" Mmax _ ^ 1) + (tl 1) 


£0 — 0 


Exemple : n=5 
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B- 2 ème cas :les niveaux équidistants : 

a. les niveaux d’énergie positive : E m — e rn > 0 
On considère le cas où E m — s m — ms 

-Dans ce cas la fonction de partition d’une seule particule s’écrit : Z a — Yïm=o e~ /}mE 


Z 


a 


n-1 


^ e -bm 
m =0 


Avec : b — (1s 

-la fonction de partition de N particules est : Z — Z C N 

-l’énergie totale : 


E 

7n 


71-1 

— V me PrnE 

z„ Z_i 


m=0 
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71-1 

M — — 'V m 2 e~ bml 

m =0 

Avec : M — — 

Ns 


-l’entropie : 

5 RE 

NK N y 


s — bM -1- Log Z a 

Avec : 

5 

S ~~NK 


Les points particuliers : (/? 

— 0, ±oo) 


o 

II 

CCI. 

-la fonction de partition 

Z a Tl 

-l’énergie totale : 

71-1 

V n — 1 

riM — } m — n - 

Lu 2 

m =0 


n — 1 

M ~ 2 

-l’entropie : 

s = logn 

Remarque : 

71 — 1 

La courbe de s=f(M) est symétrique par rapport à — j— . 


P — — 00 —> b' — —00 

-la fonction de partition 

n-1 

Z a = lim e p 2 £ 


54 





-l’énergie totale : 


-l’entropie : 


M - 


n — 1 
2 


5 — 0 


/? — +00 b — +00 


-la fonction de partition 


-l’énergie totale 


-l’entropie 


Za = 1 


M = 0 


5 — 0 


Remarque : 

La courbe de s=f(M) est symétrique par rapport à 0, car œ — 0 : 


= M b=0 - 


M-min ” 1 ” M ma x _ ^ 1 0 + (jl 1 ) 

2 ~~ 2 2 

eu = 0 


b. les niveaux d’énergie symétriques : positives pour m > 
0 et négatives pour m < 0 

L’énergie de la particule est : E m — s m — ms 

-On considère la fonction de partition d’une seule particule Z a 


w-l 

2 



71 — 1 

On fait un changement d’indice tel que : M — m + 

Donc pour m — — -* M — 0 et pour m — -> M — n — 1 
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La somme devient : 

71-1 


71-1 


•I 


n—1 \_ n-1 n \i n- l_ 1 —e n P E 

(M-—)p £ = e —(] e \ e _ M p E _ e —p E l Ç_ 


M =O 


M =O 

n-i _ l — g~ n P £ 


- e 2 


-/?£ 


2 sh((3e/ 2) 
sh(n(3e/ 2) 
shÇfie/2) 


e Ps/2 _ g-n.ps/2 


1 -e~P £ 

l — e ~npe e np£/2 _ e ~npE/2 


2sh((3s/2) 2sh(/3s/2) 


= shjnps/ 2) 
“ sh(fis/2) 


-l’énergie totale: 


E = - 


dLog Z 

dP 


= -Ai 


dLog Z a 

30 


avec Z — Z, 


N 


E = 


N dZ a 

Za dp 


sh (!j) n - c ' , ( 2 f £ ) s '*(x)' s ' , ( 2 § £ ) c ' i (x) 


Ne 

oajN 

<0 

2 



sfiiSr) 


E — — 


M 


2 

(P^ 


/nps\ (Pe\ 

n coth I —— j — coth I -—J 

zl (P^\ (npe\ 

= Ts = coth ItJ “ " cotft (—J = cot ' l(6) - " cot,l(ni,) 


M = coth(b ) — n cothÇnb) 


. .. 2E , Be 

Avec : M — — et b = — 

Ne 2 


-l’entropie 


5 = 


dKTLogZ 
dT 


= K 


^ ] [ ^£0 ] 
dT a J dB dT a J 


T —\ (—£ ) + 

KT 2 v ' 


Log z| = K[PE + Log Z] 


Or E = 


NeM b.N.M 


A 


S = NK 


S = K[ b. N. M + Log Z] = NK[b. M + Log Z a ] 


£ 

-P 2 


(npe\ (ps\ 

n coth I ■——J — coth 1 — 1 


P^ 


+ Log 


shpnpe/ 2) 


sh(Pe/2 ). 

Si N=Na (Nombre d'Avogadro) alors NK — N A K = R (constante des gaz parfaits) 
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Donc 



Pour le programme on s'intéresse a cette formule 

5 = NK[b. M + Log Z a \ 


S 

s — -= M. b + Loq Z n 

NK y a 


avec : 

5 

S -~NK 


sh(nb ) 

s = M. b + Log 

sh(b ) 


-la chaleur spécifique à volume constant : — — 

- (~ T l n coth GrO - coth (f )]) 


1 dE 
KT 2 dp 


1 Ne d 




T Ns £ 

c„ = kP t - 



+ 


1 



! f'(x') 

Avec coth (/O)) = ~ sh 2^^ = - f(x).csch 2 (f(x )) 

Et csch : cosécante hyperbolique csch(x) — s/i -1 (x) 
Donc 


= IV/é 


( t ) [ csc/ ‘ 2 ( f ) ~ " 2 csc " 2 





Si N=N a alors JV/é = = P 
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Donc 



— n 2 cosh 2 



Les points particulièrs : 


-l’entropie : 


-l’énergie totale : 


p = 0 -> b = 0 


s - logn 


M — 0 


-l’énergie totale : 


P — +00 -> b = +00 


M — +(1 — n) En utilisant le site: wolframalpha.com 


-l’entropie : 

s — 0 


Remarque : 

La courbe de s=f(M) est symétrique par rapport à 0, car m = 0 : 

.. Mmin d" M m ax _ (P- 1) + (tl 1) 

0) = M b=0 ---- 0--- 

co — 0 


C- Système à 2 niveaux application et comparaison avec l’ensemble microcanonique 


En utilisant la relation trouvée pour un système canonique à n niveau : 


Ns 

E ~ L 


(nfi£\ fps\ 

n coth I —— J — coth I -—J 


ps^ 


Pour un système à n=2 : 

On a coth(x) — —J— et x = — donc E — — — f — 7 —--J—] 

v J th(x) 2 2 Lth(2x) th(x) J 
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Et on sait que th(2x) = 


Hnrir F _ _^£ 

n ClOTiC e — 

1 +th(x) 2 2 


2(l+t/l(x) 2 ) 


-—— 1 

th(x )J 


Ns 

(1 + t/i(x) 2 ) 

1 

2 

th{x) 

th(x) 


Afe th(x) 2 

— X , 7 - = —— X th(x) 
2 th(x) 2 


-2 E 

th(x ) =- -» x = arcth 

Ns 




1 /l-hx\ 

avec :arcth(x ) = -Log ( —-J et — 1 < x < 1 


s 

X ~Y~ 2 KT 


12 K 2K 1 (Ns/2 - E\ 

-* — — —x — — x-Log[- 


1 _ K ] (Ns/2 - 

---l°g y N£ / 2 + E) 


2 

1 


\Ns/2 + Ej 
1 


et --<£<- 
2 2 


£ 1 K 

on pose s = — -* s — 2s' — — -—log 

p 2 T 2s' 6 


Finalement : 


f .Ns' - E N 
k Ns' + E y 

Ne'-E 


1 K , 

— = —r x loq 

T 2e' a Ne'+E 


c' est la même formule que pour un système micro canonique à 2 niveaux 


Donc pour N très grand, les ensembles canoniques et micro canonique donnent les mêmes 
résultats. 


IV- Progmmation en langage C : 

ce programme en C a été créer pour calculer la série de valeurs de l’entropie réduite s et 
l’énergie réduite M et enregistrer ces valeurs dans un fichier et par la suite utiliser un logiciel 
de graphisme (e.g Gnuplot) pour modéliser les résultats. 

Ce programme peut : 

-traiter le problème en utilisant les 2 ensembles; micro canonique et canonique pour n=2 et 3. 
-Pour n>3 le programme utilise l’ensemble canonique seulement, car l’ensemble micro 
canonique pose beaucoup de difficultés. 

Code source : 

g 

température négative.c 
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#include <stdio.h> 

#include <stdlib.h> 

#include <math.h> 

#include<windows.h> 

#include <conio.h> 

// IDE utiliser: Dev-C++ 5.11 

double factorial (double n); //prototype de la focntion factorial 


main ( ) 

{ 

FILE *fichier; 

f ichier=f open ( "données.txt" , "w" ); 
double M,t,b,s,Za; 

int n,ensemble=0 , répartition ,m, symétrie; 


printf ( "entrer le nombre des niveaux\nn=" ); 

scanf ( "%d" ,&n); 


//n=2 ensemble microcanonique 

if (n==2) 

{ 

printf ("\ndans le cas du système a 2 niveaux , on peut utiliser les 2 
ensembles \n 1-microcanonique (Taper 1) \n 2-canonique (Taper 2)\n"); 

scanf ( "%d" ,Sensemble); 
if (ensemble==l) 

{ 

double smin , smax, h=0 . 01 ; 
int N=1000000; 

printf ( "système a 2 niveaux de N particules(N=l.000.000)\n" ); 
fprintf(fichier,"-->M<--\t\t\t\t-->s=S/K<--\t\t\t\t-->t=KT/epsilon<-- 

\n" ) ; 

printf ("s=S/k=Nlog(N)-log((N+M)/2)*(N+M)/2-log((N-M)/2)*(N-M) /2\n") ; 

sleep (5) ; 

for (m=-N+l;m<N;m+=1500) 

// on a commencé par p=-N+l car si on commence par p=-N on va avoir un 
problème (In d'un nombre <0) dans le calcule de smin dans la lere itération 
// p=N aussi n'est pas indu , problème (In d'un nombre <0) au niveau de 
smax dans la derniere itération 

{ 

// entropie 

s= (double) N*log (N) - (double) log ( (double) (N+m) /2) * (N+m) /2- 
(double) log ( (double) (N-m) /2) * (N-m) /2 ; 

//(double) pour conversion from 'int' to 'double' , N est un entier donc 
N.log(N) va retourner un entier et pour retourner en réel on ajoute 
(double) 

printf ( "s=%f\n" , s) ; 

smax= (double) N*log(N)- (double)log ( (double) (N+m+h)/2)*(N+m+h)/2- 
(double)log ( (double) (N-m-h)/2)*(N-m-h)/2; 

smin= (double) N*log(N)- (double)log ( (double) (N+m-h)/2)*(N+m-h)/2- 
(double)log ( (double) (N-m+h)/2)*(N-m+h)/2; 

// température 

printf (" l/t=%.20fk\n" ,( smax smin) /2/h) ; —belvin 

t=pow ( (smax-smin) / (2*h) ,-1) ; c^-inverse de la pente 
printf ( "t=KT/epsilon=%.15fk\n" ,t) ; 
printf ( "%d\t%.9f\t%.9f\n" ,m, s ,t); 
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fprintf (fichier, " %d\t\t\t\t% . 91f \t\t\t\t% . 91f \n" ,m, s , t) ; 

if ( -2000<m && m<2000 ) 

//on change le pas de 1500 à 1 dans cet intervalle, car la variation de le 
température au voisinage de p=0 devient très grande 

m-=1499; 

} 

printf ( " \n\nles résultats sont enregister dans le fichier 
données.txt\n\n" ); 

} 


} 

//fin n=2 ensemble microcanonique 


//n=3 ensemble microcanonique 

if (n==3) 

{ 

printf ("\ndans le cas du système a 3 niveaux , on peut utiliser les 2 
ensembles \n\nl -microcanonique : (Taper 1) \non peut déterminer les 
contributions possibles \net calculer le nombre oméga d'états 
accessibles\n\n" ); 

printf ( "2-canonique (Taper 2)\n"); 
scanf ( "%d" , Sensemble); 
if (ensemble==l) 


{ 


FILE Contributions ; 

contributions=f open ( "contributions possibles.txt" , "w" ); 

double N,Nl,N2,N3,Nlfact,N2fact,N3fact, omega=0 , oméga m=-l, Ntfact ; 

int E ; 

fprintf (fichier, "-->M=E/e<--\t\t\t-->omega<--\t\t\t-- 
>s=S/K=log(oméga)<--\n" ); 
do{ 

printf ( "donner le nombre total des particules <181 N="); 

scanf ( "%lf" ,&N); 


} 

while (N>170); 

printf (" \nSysteme a 3 niveaux décrit par l'ensemble microcanonique 
\nles energies possibles sont: e 2e 3e\n\n"); // e c'est epsilon 

printf ("\nle nombre total des particules est: N=%.If\n\n" ,N); 
fprintf (contributions , "\nle nombre total des particules est: 

N=%.If\n\n" ,N); 

printf ("avec N= N1+N2+N3 \t et E=(N1+2N2+3N3)e\n\n" ); 

printf ("\n\n -> on va varier E de %.Ife a %.Ife \n\n" , N ,3*N); 

fprintf (contributions \n\n-> on va varier E de %.Ife a %.Ife 

\n\n" , N ,3*N); 

sleep (2); 


for (E=N;E<3*N+1;E++) 

{ 

m=omega=0; 
oméga m=-l; 

printf (" - 

- ") ; 

printf ( "\n\n\t\tl 'energie du système est: E=%de\n\n" ,E); 

fprintf (contributions - 

- \n\t\tl 'energie du système est: E=%de\n\n" ,E); 

for (N1=0;N1<N+1;N1+=1) calcul de tout les contributions 

{ 

N2=-2*N1-E+3*N; 

N3=-2*N+N1+E; 

if (N2>=0 && N3>=0) // afficher les contributions possibles 


61 













{ 

m+=l ; 

fprintf (contributions , "la contribution possible %d :\n\n",m); 

fprintf (contributions , " \t Nl=%.If N2=%.If 

N3=%.lf\n\n\n" ,N1,N2,N3) ; 

Nlfact=f actorial (NI); 

N2fact=factorial (N2); 

N3fact=f actorial (N3); 

Ntfact=Nlfact*N2fact*N3f act ; 
oméga m=f actorial (N)/Ntf act ; 
omega+=omega m; 

fprintf (contributions \tNl!=%g N2!=%g N3!=%g Nl!N2!N3=%g 

\n\n\t\tomega_%d=%g\n\n" , Nlfact ,N2f act , N3fact ,Ntf act ,m, omega_m) ; 

} 

} 


if (omega_m!=-l) // si la valeur initiale de oméga m change = 
existence des contributions possibles 

{ 

printf ( "\t==> " ); 

fprintf (contributions \t==> ") ; 

Sleep (50) ; 

printf ("le nombre des contributions possibles est: %d\n\n",m); 

fprintf (contributions , "le nombre des contributions possibles 
est: %d\n\n",m); 

printf (" \t==> "); 

fprintf (contributions \t==> ") ; 
printf ("le nombre oméga d'etats accessibe est: %g\n\n" ,oméga); 
fprintf (contributions , "le nombre oméga d'etats accessibe est: 
%g\n\n" , oméga) ; 

fprintf (fichier, "%d\t\t\t\t%g\t\t\t%.91f\n" , E , oméga , log (oméga)); 

} 

else // oméga m n'a pas changé ==> système non réalisable 

{ 

printf ("il n'existe aucune contribution qui vérifié N=%.If et 
E/e=%.d \n" , N ,E); 

fprintf (contributions , "il n'existe aucune contribution qui vérifié 
N=%.If et E/e=%.d \n" ,N,E); 

printf ("donc ce système n'est pas réalisable avec cette répartition 
des niveaux\n\n" ); 

fprintf (contributions , "donc ce système n'est pas réalisable avec 
cette répartition des niveaux\n\n" ); 

} 

} 

printf (" \n\nles résultats sont enregister dans les fichiers \n 
données.txt\n contributions possibles.txt\n\n" ); 

} 

} 

//fin n=3 ensemble microcanonique 

//ensemble canonique 

if (ensemble !=1) 

{ 

printf (" \nSysteme a %d niveaux décrit par l'ensemble canonique\n\n\n" ,n); 
printf ("ce problème sera traiter :\n" ,n); 

printf ("1er choix : des niveaux équidistants (taper l)\n",n); 
printf ("2eme choix : des niveaux non équidistants(taper 2)\n",n); 
printf ( "taper 1 ou 2\n\n",n); 

scanf ( "%d" ,Srepartition); 
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if (repartition==l) 

{ 

fprintf (fichier , "<--b=beta.e/2-->\t <--M=2E/Ne-->\t\t 

<--s=S/NK-->\n" ); // e c'est epsilon 

printf("\n\n Votre choix est équidistants (...-2e,-e,0,e,2e....)\n\n" ); 

sleep ( 1 ); 
b=n%2; 

if (b==0) 

{ 

n+=l ; 

printf("n est pair ! , pour faciliter les calculs on ajoute 1 a n pour 

le rendre pair\n la nouvelle valeur de n est n=%d\n\n" ,n); 

sleep (5) ; 

} 

printf("on va utiliser les formules suivantes\n\n" ); sleep (2); 
printf ( "M=2E/Ne=coth(b)-n.coth(nb) avec b=beta.e/2\n\n" ); sleep (3); 
printf ( "s=S/NK=Mb+log(sinh(nb)/sinh(b))\n\n" ) ; sleep (3); 

//lorsque betta=-infini 

//on considéré que la valeur -infini c'est -20 
b=-20; // <=> b --> -infiini 

s=0; //s_min 
M=n-1; //M max 

printf ( "b=%lf\tM_max=%.201f\ts_min=%.201f\n" ,b,M, s); 
fprintf (fichier , "%.3f\t\t%.201f \t\t% .201f\n" ,b,M, s); 

//calcul lorsque betta<0 

for (b=-l 8 .2+(n/10+3);b<=-0 . 2/n;b+= .01 ) 

de telle sorte que le programme peut tester les grands valeurs de n 
//après plusieurs essais on a noté que cet intervalle est la meilleure 
//car si non on va avoir soit - des nombres très grand ou bien variation 
très petite des grandeurs 

{ 

M=l/tanh(b)-n/ tanh (n*b); 

s=M*b+log(sinh(n*b)/ sinh (b)); 
printf ( "b=%lf\tM=%.201f\ts=%.201f\n" , b , M ,s); 

fprintf (fichier , "%.3f\t\t%.201f\t\t%.201f\n" ,b,M, s); 

} 

// calcul lorsque betta=0 cad b=0 

printf ( "b=0\tM=0\ts=%lf\n" , log (n)); 

fprintf (fichier, "0\t\t\t0\t\t\t%.91f\n" , log (n)); 

// calcul lorsque betta>0 

for (b=0 . 2/n;b<18 . 2-(n/10+3);b+=.01) 

{ 

M=l/tanh(b)-n/ tanh (n*b); 

s=M*b+log (sinh (n*b)/ sinh (b)); 
printf ( "b=%lf\tM=%.201f\ts=%.201f\n" ,b,M, s); 

fprintf (fichier , "%.3f \t\t% .201f \t\t% .201f\n" ,b,M, s); 

} 

//lorsque betta=+infini les details dans le fichier word programme- 
canonique . docx 
s=0; //s_min 

M=l-n; 

printf ( "b=%lf\tM_min=%.201f\ts_min=%.201f\n" ,b , M, s); 
fprintf (f ichier , .3f \t\t% .201f \t\t% .201f\n" ,b,M, s); 

} 

//fin équidistants ensemble canonique 
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else on équidistants ensemble canonique 

{ 

fprintf (fichier , "<--b=beta.e-->\t <--M=E/Ne-->\t\t <-- 

s=S/NK-->\n" ); // e c'est epsilon 

printf("\n2 possibilités pour traiter ce problème ; \nles niveau 
d'energie sont :\nl-symétriques par rapport a 0 (...-e,0,e....)\n2- 

positives (0,e,4e . ...)\n\n"); 
printf ( "taper 1 ou 2\n"); 

scanf ( "%d" ,Ssymetrie); 
if (symetrie==l) 

{ 

printf (" \nVotre choix est non équidistants et symetriques\n dans ce 
cas la répartition des niveaux d'energie suit la suite suivant\n si m>0 
Em=m A 2.e\n si non : Em=-m A 2.e\n\n" ); 

sleep (2); 
b=n%2; 

if (b==0 ) 

{ 

n+=l ; 

printf ("n est pair ! , pour faciliter les calculs on ajoute 1 a n 

pour le rendre pair\n la nouvelle valeur de n est n=%d\n\n" ,n); 

sleep (3) ; 

} 

printf ("\non va utiliser les formules suivantes\n\n" ); sleep (2); 
printf ( "M=E/Ne\n=(somme sur m de [exp(m A 2.b)]+[exp(-m A 2.b)]) A -1 * 
(somme sur m de [-m*exp(m A 2.b)]+[m*exp(-m A 2.b)]) \n avec 
b=beta.e\n\n" );sleep(4); 

printf ("s=S/NK=b.M+Log(Za) avec Za la fonction de partition d'une 
particule\n\n" ); sleep (3); 

//lorsque betta=+-infini les details dans le fichier word programme- 
canonique . docx 

//on considéré que la valeur -infini c'est -110/(n*n) 

b=(double) (60+3*n)/(n); // cette valeur est considérée comme infini 

M= (double) (n-1)*(n-1)/4 ; 
printf ("b=% .31f\tM=0\ts=0\n" ,b, -M); 
printf ("b=% .31f\tM=0\ts=0\n" ,-b,M); 

fprintf (fichier 31f\t\t0\t\t\t\t0\n" ,b,-M); 
fprintf (fichier , .31f \t\tO\ t\t\t\tO\n" , b, -M); 

// calcul lorsque betta<0 

for (b= (double) - (50+3*n) / (n) ;b<= (double) - 
0.001/ (2*n);b+= (double) 0. 01/(n) ) n<100 

{ 

Za=l;s=0; 

for (m=l;m<(n-1)/2+1;m+=l) 

{ 

Za+=exp(m*m*b) +exp (-m*m*b); 

s+=m*m*(exp(-m*m*b) -exp (m*m*b)); 

} 

M=s/Za; 

s=b*M+log(Za); 

printf ("b=%lf\tM=%.201f\ts=%.201f\n" ,b,M,s) ; 
printf ( "b=%lf\tM=-%.201f\ts=%.201f\n" , -b,M, s); 

fprintf (fi chier , "%.3f \t\t% .201f \t\t% .401f\n" ,b , M, s); 
fprintf (fichier , .3f\t\t-%.201f \t\t% .401f\n" , -b,M, s); 

} 

// calcul lorsque betta=0 cad b=0 

printf ( "b=0\tM=0\ts=%.201f\n" , (double)log (n)); 

fprintf (fichier, "0\t\t\t0\t\t\t%.91f\n" , (double)log (n)); 
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} 

// fin non équidistants et symétriques canonique 

else //non équidistants et positives canonique 

{ 

printf ("\nVotre choix est non équidistants et positives\n" ); 
printf(" dans ce cas la répartition des niveaux d'energie suit la 
suite suivant\n Em=m A 2.e\n\n" ); sleep (3); 

printf ("on va utiliser les formules suivantes\n\n" ); sleep (2); 
printf ( "M=E/Ne=(somme sur m de [exp(-m A 2.b)]) A -1 * (somme sur m de 
[m A 2.exp(-m A 2*b)]) \n avec b=beta.e\n\n");sleep(4); 

printf ( "s=S/NK=b.M+Log(Za) avec Za la fonction de partition d'une 
particule\n\n" ) ; 

sleep (3) ; 

//lorsque betta=-infini les details dans le fichier word programme- 
canonique . docx 

//on considéré que la valeur -infini c'est -110/(n*n) 

M=pow(n-1,2); 

printf ( "b=%.31f\tM=%.201f\ts=0\n" , (double) -110/(n*n),M); 
fprintf (fichier, "%.31f \t\t%lf\t\t\t\tO\n" , (double) -110/(n*n),M); 

// calcul lorsque betta<0 

for (b= (double) -(50+3*n)/(n);b<= (double) - 
0.001/(2*n);b+= (double) 0. 01/(n)) //n<100 ^ 

{ 

Za=s=0; 

for (m=0;m<n ;m+=l ) 

{ 

Za+=exp(-m*m*b); 

s+=m*m*exp(-m*m*b) ; 

} 

M=s/Za; 

s=b*M+log(Za); 

printf ( "b=%lf\tM=%.201f\ts=%.201f\n" ,b , M, s); 

fprintf (fichier , "%.3f\t\t%.201f\t\t%.401f\n" ,b,M, s); 

} 

// calcul lorsque betta=0 cad b=0 

printf ( "b=0\tM=%.201f\ts=%.201f\n" , (double) 1/6*(n-1)*(2*n- 
1), (double)log (n)); 

fprintf (fichier, "0\t\t\t%.91f \t\t\t% .91f\n" , (double) 1/6*(n- 
1)* (2*n-l), (double)log (n)); 

// calcul lorsque betta>0 

for (b= (double) 0.001/(2*n);b<= (double) (50+3*n)/(n);b+= (double) 0.01/(n) 

) //n<100 

{ 

Za=s=0; 

for (m=0;m<n;m+=l) 

{ 

Za+=exp(-m*m*b) ; 

s+=m*m*exp(-m*m*b); 

} 

M=s/Za; 

s=b*M+log(Za); 

printf ( "b=%lf\tM=%.201f\ts=%.201f\n" ,b , M ,s); 

fprintf (fi chier , "%.3f \t\t% .201f \t\t% .401f\n" ,b,M, s); 

} 
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//lorsque betta=+infini les details dans le fichier word programme- 
canonique . docx 

printf ( "b=%.31f\tM=0\ts=0\n" ,b= (double) (60+3*n)/(n)); 
fprintf (fichier, "% .31f \t\t%lf\t\t\t\tO\n" ,b= (double) (60+3*n)/(n)) ; 

} 

//fin non équidistants ensemble canonique 

} 

//fin non équidistants choix de la symétrie ensemble canonique 

printf (" \n\nles résultats sont enregister dans le fichier 
données.txt\n\n" ); 

} 

// fin ensemble canonique 

fclose (fichier); 

System ("PAUSE") ; 

} 

//focntion factorial 

double factorial (double n) 

{ 

if (n>l) 

return n*f actorial (n-1); 

else 

return 1 ; 

} 


66 







CHAPITRE 6 : Champ d’application 


Maintenir de la matière à une température absolue négative peut aboutir à des applications 
étonnantes. 


l.en thermodyamique : 

es températures négatives, permettent en théorie des moteurs thermique selon le cycle de 
Carnot, dont le rendement dépasse 100%. 

On caractérise la performance des systèmes moteurs par leurs 
rendements R. 

Ce rendement est définie comme le rapport du travail produit sur 
la quantité de 
Chaleur consommée, soit : 


q = e = 7 


énergie utile 


énergie fournie machine 


W 

0 ^ 


w 

* 


Ces systèmes fonctionnent selon un cycle de CARNOT 
Le cycle de Carnot est par définition réversible: 


W+Qc+Qf =0 


Q£ + Qf =0 

Te Tf 


% _ -W _ Qc+Qf _ , , Qf 

Cmoteui— — Z — i + ■ 


Qc 


Qc 


Qc 


Qc 

Te 


v.. 

Tf Qc 


IL 

Te 



Représentation du cycle : 



Qf t f 

D OU ecamot— 1 + ~~ = 1 

Qc Te 

Le rendement e ca mot est toujours <1 pout les systèmes à température positive. 

Par contre si Tf et Tc ne sont pas de même signe (Tf < 0, T c > 0), le rendement devient 
supérieur à l'unité. [ - ] 


2.en cosmologie : 

Lorsqu’un un gaz d’atomes ultra-froids stabilisé à quelques fractions de kelvin en dessous du 
zéro absolu, il devrait normalement s’effondrer sur lui-même (similaire à 
une pression négative). En effet, le passage à une température absolue négative s’accompagne 
d’un changement de signe dans de nombreuses quantités comme l’interaction entre les 
particules. Ces particules qui se repoussaient auparavant, s'attirent désormais (de la même 
façon qu’une étoile s’effondre sur elle-même sous l’effet de la gravitation à la fin de sa vie car 
les réactions thermonucléaires qui la stabilisaient s’arrêtent), 

Rien n'est supposé pouvoir empêcher l'effondrement du gaz! En fait, c'est la pression qui vient 
nous sauver la mise! On peut montrer avec très peu de calculs que ce système à température 
absolue négative n'est stable que si la pression devient négative également. Ces effets 
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'absurdes' observés dans l'expérience des chercheurs allemands sont en effet par certains 
aspects similaires à ce qui est observé dans le cadre de l'énergie noire. 

Donc la température négative pourrait aussi être utile en cosmologie, car le comportement de 
la matière à une température absolue négative présente des parallèles avec la mystérieuse 
énergie noire qui permet l’expansion de l’univers. 


Les pressions négatives 

Afin d'étudier les conditions d'équilibre thermique, on considère un gaz dans une boîte du 
volume Vbox, avec une énergie totale E fixe. Dans l'équilibre thermique l'entropie S prend la 
valeur maximale possible sous les contraintes données, on note que le volume de la boîte est 
seulement une limite supérieure pour le volume V du gaz (V < Vbox). 

Dans la plupart des cas, ces conditions conduisent à des pressions positives P> 0, ce qui 
oblige le gaz à remplir toute la boîte (V = Vbox). 

De façon plus générale, cependant, le principe d'entropie maximale ne nécessite que: 

^l E > 0 

dV 

si cette dérivée est négative le système peut se contracter spontanément, et ainsi, augmenter 
son entropie, en conséquence, il n'y aurait pas de solution à l'équilibre avec V=Vb ox , mais le 
système serait instable contre l'effondrement. 

D'après la différentielle totale de l'energie interne: 

dE=TdS - pdV => ds = - dE- - dv 

r TT 

donc: |^Ie = - or ^-Ie > 0 donc - > 0 
dv T dv T 

V P 

A l'équilibre il faut que la pression absolue et la température soient de même signe - > 0. 
Donc les températures négatives impliquent les pressions négatives et vice versa. 
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CONCLUSION 


A partir des donnés de la théorie de Norman Ramsey présentée en 1956 et de l’expérience des 
allemands concernant les moteurs thermiques dont le rendement est supérieur à l’unité, on 
peut tirer que ceci est valable seulement pour les systèmes qui ont déjà des températures 
absolues négatives, car si c’est nous qui les produits, on doit toujours fournir de l’énergie au 
système, pour ne pas le laisser refroidir, donc l’alimenter en permanence, et donc la sur unité 
gagnée du moteur est largement faible par rapport à l’énergie qu’il faut fournir afin que ce 
système reste à température négative . 

Donc, la seule utilité vraie de ce genre de système serait que ce système soit déjà tout seul 
dans ce mode de température, sans qu'on ait besoin de l'alimenter pour cela, sans lui fournir de 
l'énergie. C'est le cas de l'énergie noire si on la considère comme un gaz d'énergie partout 
présent, dans ce cas on ne fera que dire qu'on peut récupérer de l'énergie depuis une source 
existant dans l'univers. 

Tout ça reste une revue scientifique pour faire une révolution scientifique qui bouleverse le 
monde entier, car jusqu’à maintenant on ne sait pas comment détecter cette énergie noire. 
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